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The  concerns  of  energy  consumption  and  environment  pollution  urge  researchers  to  work  on  the 
development  of  clean  energy  and  the  utilization  of  waste  energy.  As  one  important  topic,  absorption 
cycle  technology  has  attracted  considerable  attention  because  it  can  be  powered  by  low-grade  heat,  e.g., 
solar  energy  and  waste  heat.  In  recent  years,  researchers  proposed  that  ionic  liquids  (ILs)  as  novel 
alternative  absorbent  combined  with  refrigerant  such  as  water,  ammonia,  alcohols,  and  hydrofluor¬ 
ocarbons  can  be  used  as  working  pairs  for  absorption  refrigeration  cycle,  heat  pump,  and  absorption 
power  cycle.  In  this  paper,  researches  done  in  imidazolium  IL  working  pairs  regarding  to  status  of 
evaluation  and  selection  methods,  thermophysical  property  measurement  and  modeling,  as  well  as  their 
future  prospect  assessments,  i.e.,  developing  potential  studies  about  the  absorption  cycle  performance 
adopting  new  working  pairs,  have  been  reviewed. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

2.2.  Significance  ofR&D  of  absorption  cycle  and  ionic  liquid  working 
pair 

To  deal  with  the  changes  in  global  climate,  rapid  growth  of 
energy  demand  and  rapid  development  of  renewable  energy  and 
energy  utilization,  much  attention  has  been  focused  on  the 
research  and  development  (R&D)  of  absorption  cycles,  such  as 
absorption  heat  pump  and  absorption  refrigeration  cycle,  because 
they  can  make  effective  use  of  low-grade  heat  [1-4],  However,  the 
R&D  of  absorption  technology  faces  a  series  of  problems.  For 
example,  the  temperature  range  of  the  heat  source  is  wider,  with 
a  range  of  80-90  C  collected  by  economic  and  efficient  flat  solar 
collectors  to  500  -  600  °C  exhausted  from  gas  turbines  or  internal 
combustion  engines.  The  functions  and  configurations  of  new 
absorption  cycles  are  more  diverse  and  complex,  such  as  cogen¬ 
eration  systems  that  combine  power  and  refrigeration  cycles  [5-7], 
compression  and  absorption  hybrid  cycles  driven  by  low-grade 
heat  and  electrical  energy  [8,9],  and  distributed  energy  systems 
that  combine  cooling,  heating,  and  power  cycles  [10]  involving 
energy  storage  [11-15]  and  absorption  technologies  [16-18],  The 
above  developments  present  new  demands  for  the  R&D  of  work¬ 
ing  pairs  of  absorption  cycles  as  a  key  and  essential  work. 

Generally,  the  working  pair  of  absorption  cycle  is  a  multi¬ 
species  solution,  which  is  usually  a  binary  solution  consisting  of 
two  species  with  a  large  boiling  point  difference.  The  species  with 
relatively  low  boiling  point,  i.e.,  the  volatile  species,  acts  as 
refrigerant,  whereas  the  species  with  high  boiling  point,  i.e.,  the 
non-volatile  species,  acts  as  absorbent.  The  main  types,  existing 
problems,  and  requests  for  future  development  of  the  absorption 
cycle  working  pairs  are  briefly  summarized  and  listed  in  Table  1. 
Based  on  the  species  of  refrigerant,  the  working  pairs  are  divided 
into  the  following  five  categories:  water,  ammonia,  alcohol, 
halohydrocarbon  and  hydrocarbon  (HC).  The  combinations  of  five 
kinds  of  refrigerants  and  different  absorbents  can  basically  meet 
the  needs  of  different  applications.  For  example,  the  H20/LiBr 
working  pair  is  mainly  applied  in  room  air  conditioning.  The  NH3/ 
H20  working  pair  can  meet  cooling  requirements  below  0  C.  The 
two  working  pairs  described  above  are  the  most  widely  used 
systems  with  environmentally  friendly  features. 


However,  the  H20/LiBr  system  presents  certain  disadvantages, 
such  as  crystallization,  corrosion,  and  negative  pressure  operation. 
Some  researchers  have  attempted  to  improve  the  H20/LiBr  work¬ 
ing  pair  by  adding  auxiliary  species  [19,20],  The  NH3/H20  system 
also  has  some  drawbacks,  such  as  difficulty  in  separation  and 
toxicity.  Some  researchers  have  added  non-volatile  species  to 
reduce  the  separation  difficulty  [21-23], 

Alcohol  working  pairs  are  mainly  composed  of  low  carbon 
alcohols,  such  as  methanol  and  ethanol,  or  fluoroalcohols,  such  as 
2,2,2-trifluoroethanol  (TFE)  and  hexafluoroisopropanol  (HFIP),  as 
refrigerant  and  salts  or  high-boiling  point  organics  as  absorbent. 
Alcohols  are  inflammable  and  methanol  is  toxic,  so  greater  interest 
has  been  given  to  fluoroalcohols.  Fluoroalcohols  are  non-corrosive, 
non-combustible,  and  have  good  thermal  characteristics.  TFE  com¬ 
posed  of  working  pairs  with  some  high-boiling  point  organics  have 
been  used,  such  as  tetraethylene  glycol  dimethyl  ether  (E181 )  [24,25], 
N-methyl-2-pyrrolidone  (NMP)  [26],  N,N'-dimethylpropyleneura 
(DMPU)  [27],  N,N'-dimethylethyleneurea  (DMEU)  [27],  and  quinoline 
[28,29],  These  working  pairs  generally  have  no  crystal  restrictions, 
and  therefore,  have  relatively  wide  operating  ranges. 

Hydrofluorocarbons  (HFCs)  are  refrigerants  with  benign  tech¬ 
nical  properties,  which  are  extensively  used  in  the  vapor  compres¬ 
sion  refrigeration  cycle.  It  is  worth  mentioning  that  HFCs  used 
in  the  absorption  refrigeration  cycle  are  almost  harmless  to  the 
atmospheric  ozone,  whereas  chlorofiuorocarbons  (CFCs)  are 
harmful  to  the  atmospheric  ozone  commonly  cannot  be  used  for 
the  absorption  refrigeration  cycle  [30].  HC  working  pairs  have  also 
attracted  widespread  attentions  because  of  their  naturally  eco- 
friendly  feature  [8,31],  However,  due  to  the  flammability  and 
explosibility,  researches  and  applications  of  HC  working  pairs  are 
impeded  by  numerous  restrictions. 

Both  H20/LiBr  and  NH3/H20  systems  can  be  developed  for 
future  applications,  but  existing  problems  need  to  be  solved 
urgently.  Unlike  H20  and  NH3,  absorbents  suitable  for  alcohols, 
HFCs,  and  HCs  are  still  quite  few.  Considering  that  HFCs  and  HCs 
can  meet  the  refrigeration  technology  demands  of  H20  and  NH3, 
absorption  cycle  working  pairs  containing  HFC  or  HC  have  great 
potential  for  exploration  and  are  expected  to  fill  the  gaps  of  H20/ 
LiBr  or  NH3/H20  systems.  Around  exploring  advanced  absorbents 
with  HFC  and  HC,  researchers  are  exerting  efforts  to  develop  novel 
absorption  cycle  working  pairs  [27,32], 


Table  1 

R&D  of  absorption  cycle  working  pairs. 


Refrigerant 

species 


Absorbent 

species 


Existing  problems 


Main  of  the  past  R&D  Requests  of  future  development 


h2o 

LiBr 

Corrosion,  crystallization 

Adding  other  auxiliary  species,  e.g., 
octanol  and  other  species 

Abate  or  avoid  the  problems  of  corrosion  and 
crystallization 

nh3 

H20 

Toxicity,  difficulty  in  separating  NH3  from  H20 

Other  organic  ammoniates 

Reduce  the  separating  energy  consumption 

1 )  Absences  of  techno-economic 

1 )  Meet  needs  of  the  R&D  of  new  refrigerants, 

HFC  &  HC 

Organic 

solvents 

significance 

2)  Values  of  ODPa  and  GWPb  cannot  meet 
the  needs  of  developing 

New  refrigerants  with  higher  ODP 
and  GWP  values 

e.g.,  R152a,  R32 

2)  Suitable  for  novel  hybrid  cycles  and  co¬ 
generation  cycles 

Alcohol  (CH3OH, 

Organic 

To  be  devoid  of  techno-economic  significance 

Few 

Better  techno-economic  significance,  and 

C2H5OH,  TFE) 

solvents 

as  like  as  the  system  H20/LiBr 

environmentally  friendly  characteristics 

Ozone  depletion  potential. 
Global  warming  potential. 
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3.2.  Basic  properties  and  characteristics  of  ILs  used  as  absorbent  (5)  The  chemical  and  physical  properties  of  ILs  can  be  adjusted  by 

species  the  design  of  anion  and  cation. 


ILs  are  composed  of  ions,  and  generally  consist  of  organic 
cations  containing  nitrogen  or  phosphorus  and  organic  or  inor¬ 
ganic  anions  [33],  In  1992,  Wilkes  et  al.  [34]  published  reports  on 
synthesis  and  performance  of  l-ethyl-3-methylimidazolium  tetra- 
fluoroborate  ([EMIM]BF4).  They  found  that  this  IL  is  moisture 
stable,  does  not  undergo  hydrolysis,  has  a  low  melting  point,  and 
has  good  thermal  stability.  Thereafter,  the  study  of  ILs  has  been 
greatly  promoted,  constantly  developed,  and  emerged  worldwide 
upsurge. 

Variety  species  of  ILs  exist,  and  their  physicochemical  proper¬ 
ties  can  be  regulated  by  changing  the  type  and  structure  of  anions 
and  cations.  Thus,  ILs  are  called  designable  solvents  [35,36], 
Generally,  based  on  the  types  of  organic  cation,  ILs  can  be  divided 
into  imidazolium,  pyridinium,  tetraalkyl  ammonium,  tetraalkyl 
phosphonium,  and  so  on  [37,38],  Common  cationic  structures 
are  shown  in  Fig.  1.  Among  ILs  mentioned  above,  imidazolium 
ILs  are  the  most  stable,  most  common,  and  most  studied. 

According  to  the  solubility  differences  of  ILs  and  water,  ILs 
can  be  divided  into  two  categories:  hydrophilic  ILs,  such  as 
l-alkyl-3-methylimidazolium  chloride  ([CnMIl\/I]Cl,  n  is  the 
number  of  C  atom  in  the  alkyl,  n  =  l,  2,  3,...),  l-alkyl-3-methyl- 
imidazolium  tetrafluoroborate  ([CnMIM]BF4),  l-alkyl-3- 
methylimidazolium  dimethylphosphate  ([C„MIM]DMP),  1- 
alkyl-3-methylimidazolium  methylsulfate  ([CnMIM]MeS04), 
and  l-alkyl-3-methylimidazolium  ethylsulfate  ([CnMIM]Et 
S04),  and  hydrophobic  ILs,  such  as  l-alkyl-3-methylimidazo- 
lium  hexafluorophosphate  ([CnMIMJPF6),  l-alkyl-3-methyl- 
imidazolium  trifluoromethanesulfonate  ([CnMIM]TfO),  and  1- 
alkyl-3-methylimidazolium  bis(trifluoromethylsulfonyl)imide 
([CnMIM]Tf2N)  [39,40].  The  hydrophilicity  and  hydrophobicity 
of  ILs  are  affected  by  type  and  structure  of  anions  and  cations. 
Usually,  a  shorter  alkyl  side  chain  leads  to  stronger  hydrophi¬ 
licity  of  IL,  conversely,  a  longer  alkyl  side  chain  leads  to 
stronger  hydrophobicity  of  IL. 

As  a  new  kind  of  environmentally  friendly  solvents,  ILs  have 
been  widely  utilized  in  chemical  reactions,  catalysis,  separation, 
and  electrochemistry  [36,40,41],  Moreover,  ILs  can  be  used  as 
absorbent  species  in  combination  with  refrigerants  to  constitute 
new  working  pairs  for  absorption  cycle.  Compared  with  organic 
solvents,  the  advantages  of  IL  are  mainly  embodied  in  the  follow¬ 
ing  aspects. 

(1)  Most  ILs  exist  as  liquid  at  a  wide  range  around  room  tem¬ 
perature.  The  boiling  points  of  ILs  are  much  higher  than  those 
of  refrigerants  because  of  the  ignorable  vapor  pressure.  In  the 
generation  process,  the  refrigerant  can  be  easily  separated 
from  IL  with  high  purity  in  the  cycle. 

(2)  The  heat  capacity  of  IL  solution  is  small,  which  is  beneficial  in 
improving  the  cycle  efficiency. 

(3)  Solubilities  with  inorganic  or  organic  species  and  the  affinity 
with  refrigerants  are  good,  which  is  in  favor  of  enhancing  the 
mass  transfer. 

(4)  ILs  have  high  chemical  and  thermal  stability,  high  thermal 
decomposition  temperature,  and  are  non-flammable. 


R4  R, 

R,  N  *  R.  p- 

R/  R3  R/  R3 


Imidazolium  Pyridinium  Tetraalkyl 

ammonium  phosphomum 

Fig.  1.  Structure  of  the  cations  of  the  common  ILs. 


Many  research  groups  have  investigated  the  thermophysical 
properties  of  working  pairs  consisting  of  ILs  and  common  refrig¬ 
erants,  such  as  H20,  NH3i  alcohol,  and  HFCs,  for  absorption  cycle 
[42-45].  These  researches  exhibit  the  great  potential  and  signifi¬ 
cant  application  prospects  of  these  novel  working  pairs. 

3.3.  Previous  studies 

In  2004,  Kim  et  al.  [42]  proposed  some  working  pairs  composed  of 
imidazolium  ILs  (l-butyl-3-methylimidazolium  tetrafluoroborate 
([BMIM]BF4),  l-butyl-3-methylimidazolium  bromide  ([BMIMJBr)) 
and  refrigerants  (TFE,  1,1,1,2-tetrafluoroethane  (R134a)).  Among  these 
systems,  TFE/[BMIM]Br  and  TFE/[BMIM]BF4  were  selected  as  alter¬ 
native  working  pairs  for  absorption  refrigeration  cycle  and  heat 
pumps.  The  investigated  results  of  heat  capacity  and  vapor  pressure 
of  the  two  systems  show  that  [BMIM]Br  and  TFE  had  a  stronger 
interaction,  and  IL  working  pairs  generally  have  good  thermodynamic 
properties  and  deserve  further  study. 

Since  2006,  Shiflett  and  Yokozeki  [46-48]  from  DuPont  Com¬ 
pany  have  published  lots  of  patents  and  papers  on  the  use  of  ILs  as 
absorbents,  which  are  impressive  researches  at  this  field  in  recent 
decade.  The  solubilities  and  diffusivities  of  H20  [49],  NH3  [23,50], 
C02  [44,51]  and  HFCs  [52-62]  in  ILs  were  measured,  and  some 
new  working  pairs  for  absorption  cycle  were  proposed,  such  as 
H20/[EMIM]BF4,  NH3/[DMEA]Ac  (N,N-dimethylethanolammonium 
acetate),  and  R134a/[EMIM]BEI  (l-ethyl-3-methylimidazolium  bis 
(pentafluoroethylsulfonyl)imide). 

In  2006,  Sen  and  Paolucci  [45]  presented  that  IL  l-butyl-3- 
methylimidazolium  hexafluorophosphate  ([BMIM]PF6)  can  be 
used  in  absorption  refrigeration  cycle  with  C02.  Some  properties, 
such  as  vapor  pressure,  thermal  stability,  and  dissolve  ability,  were 
discussed,  which  demonstrated  that  thermodynamic  data  of  IL 
working  pair  systems  are  still  lacking.  In  2010,  Martin  et  al.  [63] 
selected  appropriate  IL  and  supercritical  C02  as  working  pairs  and 
calculated  the  coefficient  of  performance  (COP)  for  absorption 
refrigeration  cycle  based  on  the  phase  equilibrium  of  the  system. 
However,  because  of  the  necessity  of  operating  with  a  higher 
circulation  ratio,  the  COP  was  lower  than  that  of  the  conventional 
system  NH3/H20.  In  2012,  Kim  et  al.  [64,65]  conducted  theoretical 
researches  about  the  thermodynamic  performance  of  a  miniature 
absorption  refrigeration  system  with  refrigerant/IL  as  working 
pairs.  Previously  studied  ILs  including  l-ethyl-3-methyli- 
midazolium  bis(trifluoromethylsulfonyl)imide  ([EMIM]Tf2N), 
[EMIM]BF4,  [BMIM]BF4,  [BMIM]PF6,  l-hexyl-3-methylimidazo- 
lium  bis(trifluoromethylsulfonyl)imide  ([HMIM]Tf2N),  l-hexyl-3- 
methylimidazolium  tetrafluoroborate  ([HMIM]BF4),  and  1-hexyl- 
3-methylimidazolium  hexafluorophosphate  ([HMIM]PF6),  as  well 
as  refrigerants  including  H20,  pentafluoroethane  (R125),  R134a, 
1,1,1 -trifluoroethane  (R143a),  1,1-difluoroethane  (R152a),  difluor- 
omethane  (R32),  l,2-dichloro-l,l,2,2-tetrafluoroethane  (R114),  and 
1,1,2,2-tetrafluoroethane  (R134).  COP  and  feasibility  were  calcu¬ 
lated  and  compared.  Results  showed  that  among  these  systems, 
H20/[EMIM]BF4  had  the  highest  COP  with  a  value  of  0.91. 

Liang  et  al.  [66]  and  Zhao  et  al.  [67]  determined  the  heat 
capacity,  viscosity,  density  and  phase  equilibrium  for  methanol/ 
[BMIMjCl  (l-butyl-3-methylimidazolium  chloride)  and  methanol/ 
[DM1M]DMP  (1,3-dimethylimidazolium  dimethylphosphate)  sys¬ 
tems.  They  also  analyzed  the  potential  of  the  two  systems,  which 
were  found  to  meet  the  requirements  of  benign  working  pairs  of 
absorption  cycle. 

In  2011,  Zhang  et  al.  [68]  evaluated  the  performance  of  the  H20/ 
[EMIMJDMP  (l-ethyl-3-methylimidazolium  dimethylphosphate) 
system  in  absorption  refrigeration  cycle,  results  indicated  that  its 
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COP  was  lower  than  that  of  the  traditional  H20/LiBr  system  but 
still  higher  than  0.7.  Moreover,  the  generation  temperature  of  the 
H20/[EM1M]DMP  system  was  lower  than  that  of  the  H20/LiBr 
system,  showing  that  the  new  cycle  could  be  driven  by  waste  heat 
or  hot  water  with  lower  temperature.  Zhang  et  al.  [69,70]  also 
simulated  the  thermodynamic  performance  of  the  heat  pump 
using  H20/[EMIM]DMP  and  H20/[DMIM]DMP  systems  as  working 
pairs,  and  compared  the  result  with  that  of  H20/LiBr  and  TFE/E181 
systems.  The  results  indicated  that  the  COP  of  H20/[EMIM]DMP 
and  H20/[DM1M]DMP  systems  were  slightly  lower  than  that  of  the 
H20/LiBr  system,  but  significantly  higher  than  the  value  of  the  TFE / 
E181  system. 

Since  2007,  author's  group  has  conducted  numerous  theoretical 
and  experimental  studies  on  1L  working  pairs,  and  published 
several  related  papers  [71-83],  Based  on  the  solution  thermody¬ 
namic  theory,  the  excess  Gibbs  function  (CE)  evaluation  method  is 
expanded  to  IL  working  pairs  using  the  UNIFAC  model  as  calcula¬ 
tion  tool,  and  this  method  is  applied  in  the  screening  of  H20/IL, 
NH3/IL,  and  HFC/IL  systems.  According  to  the  evaluation  method, 
hydrophilic  IL  1,3-dimethylimidazolium  chloride  ([DM1M]C1)  was 
selected  and  synthesized.  A  variety  of  experimental  methods  were 
adopted  to  determine  thermophysical  properties,  such  as  density, 
heat  capacity,  vapor  pressure  and  solubility,  for  H20/IL,  NH3/IL,  and 
HFC/IL  systems.  Appropriate  thermodynamic  models  were  chosen 
to  regress  and  predict  the  properties  of  selected  systems.  The 
influence  of  ILs  as  additives  in  traditional  H20/LiBr  and  NH3/H20 
systems  was  investigated.  Furthermore,  the  characteristics  of 
absorption  refrigeration  cycle  and  heat  pump  using  IL  working 
pairs  were  simulated  and  calculated. 

Refrigerant/IL  systems  have  increasingly  received  much  atten¬ 
tion  as  novel,  environmentally  friendly  working  pairs.  Researches 
show  that  IL  working  pairs  can  extend  the  operating  temperature 
range  of  the  cycle  and  strengthen  the  separation  and  generation 
processes  of  absorbate,  thus  reducing  the  energy  consumption  and 
potentially  reducing  or  eliminating  corrosion  and  crystallization. 
Proper  combinations  of  refrigerant  and  IL  may  generate  new 
working  pairs  with  excellent  performance.  However,  some  pro¬ 
blems  with  IL  working  pairs  are  still  encountered,  such  as  high 
circulation  ratio  and  low  energy  conversion  efficiency,  showing 
the  necessity  for  further  development.  Fortunately,  ILs  have  huge 
species  number,  and  their  functions  can  be  designed  by  changing 
the  structure  of  cation  and  the  species  of  anion  to  meet  the 
increasing  demands  for  new  working  pairs.  Therefore,  selecting 
more  valuable  working  pairs  from  a  large  number  of  research 
objects,  measuring  necessary  chemical  and  physical  properties  of 
new  systems  with  high  accuracy,  establishing  corresponding 
models  of  physical  and  chemical  properties,  and  evaluating  the 
application  significances  of  new  IL  working  pairs,  are  still  challen¬ 
ging  issues  that  must  be  overcome  to  accelerate  the  application 
progress  of  IL  working  pairs. 

3.4.  Focus  and  proposals  of  this  work 

This  work  aims  to  exhibit  the  international  research  frontier  of 
absorption  cycle  IL  working  pairs  in  recent  years  to  make  relevant 
researches  more  rational  and  efficient.  The  evaluation  method  for 
IL  working  pairs  was  introduced,  which  was  based  on  macroscopic 
properties  and  intermolecular  interactions,  combining  the  UNIFAC 
model  with  extreme  CE  criterion.  The  research  progress  on  the 
thermophysical  properties  of  IL  working  pairs  was  presented, 
including  the  measurement  method  and  modeling  of  vapor 
pressure,  solubility,  heat  capacity  and  density  for  popular  systems. 
Particularly,  on  the  basis  of  researches  on  thermophysical  proper¬ 
ties  of  systems  containing  refrigerants  H20,  NH3,  and  HFCs,  several 
binary  and  ternary  systems  with  potential  have  been  presented. 
Finally,  the  performance  of  absorption  cycle  using  IL  working  pairs 


was  calculated  and  evaluated  to  verify  further  the  future  develop¬ 
mental  potential  of  new  systems. 


2.  Selection  of  IL  absorbent  for  working  pair  innovation 

2.1.  General  behaviors  of  the  vapor-liquid  equilibrium 

The  activity  coefficient  of  the  refrigerant  in  IL  can  be  calculated 
using  the  vapor-liquid  equilibrium  (VLE)  data  of  IL  working  pairs 
and  can  be  used  to  quantitatively  characterize  the  volatility  of 
refrigerant  and  the  solubility  in  IL,  establish  the  thermodynamic 
model  of  IL  working  pair,  analyze  the  influence  of  component 
structure  on  the  interaction  between  refrigerant  and  IL,  and  guide 
the  selection  and  optimization  of  IL  absorbent. 

For  refrigerant  (1)/IL  (2)  binary  systems,  the  concentration  of  IL 
in  the  gas  phase  can  be  neglected,  i.e.,  the  gas  phase  mole  fraction 
of  the  refrigerant  species  is  equal  to  1.  Hence,  the  VLE  relationship 
can  be  expressed  as  [84] 

p(T)  =  x,/1ps1(T)  (1) 

where  p  is  the  system  pressure,  Xi,  }q,  and  p\  are  the  liquid  phase 
mole  fraction,  activity  coefficient,  and  saturated  vapor  pressure  of 
the  refrigerant  species,  respectively.  Generally,  the  activity  coeffi¬ 
cient  model,  Wilson  model,  UNIQUAC  model,  NRTL  model,  and  so 
on  can  be  adopted  to  correlate  the  jq  based  on  the  experimental 
VLE  data  of  working  pairs.  Appropriate  model  is  selected  by 
considering  several  aspects  of  theoretical  analysis  and  experimen¬ 
tal  measurement.  Therefore,  the  activity  coefficient  value  of  the 
refrigerant  species  jq  can  be  correlated  and  predicted,  while  the 
activity  coefficient  value  of  the  IL  species  y2  cannot  be  calculated 
directly,  because  ILs  have  negligible  vapor  pressure.  However, 
when  jq  can  be  estimated,  y2  can  be  determined  by  integration 
of  the  Gibbs-Duhem  equation  at  constant  temperature  and  pres¬ 
sure  [85], 

din  y-l/dx2+x2  din  y2/dx2  =  0  (2) 

where  x2  is  the  liquid  phase  mole  fraction  of  the  IL  species.  For 
instance,  the  activity  coefficient  of  refrigerant  species  can  be 
represented  by  an  appropriate  model  [85] 

In  Xi  =  ox2  +  bx2  +  cx2  (3) 

where  a,  b,  and  c  are  concentration-independent  parameters.  The 
value  of  In  y2  for  the  IL  species  can  then  be  obtained  in  terms  of 
the  same  parameters. 

In  addition,  the  partial  pressure  of  Raoult's  law  used  for 
describing  the  VLE  behavior  of  the  refrigerant  in  a  liquid  ideal 
mixture  is  expressed  as  [84] 

PiL(7")  =  XiPj(T)  (i=l,2,  ...,N)  (4) 

In  a  p  -  x(y )  diagram,  Eq.  (4)  exhibits  a  straight  line.  The  p  -  x(y ) 
curve  of  the  non-ideal  mixtures  that  display  positive  deviation 
from  Raoult’s  law  falls  above  the  line;  whereas  the  p-x(y)  curve  of 
the  non-ideal  mixtures  that  show  negative  deviation  from  Raoult's 
law  falls  below  the  line,  causing  the  negative  deviation  of  the 
liquid  phase  from  the  ideal  mixtures,  exhibiting  the  value  of  yt  is 
less  than  1. 

In  the  view  of  vapor  pressure  criterion  for  selecting  absorption 
cycle  working  pairs,  systems  with  strong  absorbing  ability  are 
usually  classified  as  exhibiting  negative  deviation  from  Raoult's 
law  [86],  Working  pairs  exhibiting  highly  negative  deviation  from 
Raoult's  law  produce  the  best  results,  as  Morrissey  has  described 
[87],  The  solvent  vapor  pressure  decreases  because  the  solvent 
molecules  in  the  liquid  phase  have  strong  affinity  to  the  solute 
molecules. 
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Fig.  2.  Vapor  pressure-composition  diagram  of  several  IL  working  pairs.  ■,  H20 / 
[DM1MJDMP,  323.2  K  [72];  \  TFE/[EM1M]BF„,  323.15  I<  [80];  •,  R32/]EMIM]Tf2N, 
283.15  I<  [54];  t,  NH3/]BMIM]PF6,  283.4  K  [23],  Dashed  line:  Raoult's  law. 


Table  2 

Bond  energy  and  bond  length  of  some  hydrogen  bonds  [86], 


Hydrogen  bond 

Bond  energy/kj  mol  1 

Bond  length/pm 

Compounds 

F-H 

•F 

28.1 

255 

(HF)„ 

0— H 

0 

18.8 

276 

H20  (ice) 

0— H 

0 

25.9 

266 

CH3OH,  C2H5OH 

N-H-  • 

•F 

20.9 

268 

nh4f 

N-H-  • 

0 

20.9 

286 

ch3conh2 

N-H  -  • 

•N 

5.4 

338 

nh3 

Studies  on  available  IL  working  pairs  of  binary  and  ternary 
systems  and  common  refrigerants  H20,  NH3,  alcohols,  HCs,  and 
HFCs  have  shown  that  the  VLE  behaviors  of  IL  working  pairs 
generally  manifest  as  non-ideal  mixtures.  For  example,  the  vapor 
pressure-composition  (p-x)  diagram  in  Fig.  2  describes  the 
behavior  of  several  IL  working  pairs  showing  negative  deviation 
from  Raoult’s  Law. 

2.2.  Effects  of  the  molecular  structure 

Essentially,  the  non-ideal  VLE  behavior  of  the  working  pair 
depends  on  intermolecular  force,  which  includes  the  interaction 
between  different  molecules  and  homologous  molecules  for  a 
multi-species  system.  The  quantitative  relationships  for  contacting 
intermolecular  forces  and  macroscopic  properties  are  now  con¬ 
fined  to  simple  idealized  systems.  The  mechanism  underlying 
solubility  difference  among  IL  working  pairs  is  not  clear  at  the 
molecular  level,  the  hydrogen  bond  is  only  known  to  play  an 
important  role  [62,88], 

The  hydrogen  bond  can  be  simply  expressed  as  X— H  ■  ■■  Y, 
where  the  dotted  line  represents  the  hydrogen  bond.  Based  on  the 
formation  theory  of  hydrogen  bond,  the  strength  of  the  hydrogen 
bond  is  associated  with  the  electronegativity  of  X  and  Y  atoms 
connected  with  the  H  atom.  The  stronger  the  electronegativity  of 
X  and  Y  atoms  is,  the  stronger  the  hydrogen  bond  will  be. 
Meanwhile,  the  hydrogen  bond  is  also  influenced  by  the  radius 
of  the  Y  atom.  A  smaller  the  radius  of  the  Y  atom  results  in  a 
stronger  hydrogen  bond.  The  bond  energy  and  bond  length  of 
some  hydrogen  bonds  are  listed  in  Table  2  [89],  Some  atomic 
electronegativity  values  (Pauling  scale)  are  shown  in  Fig.  3  [90].  In 
the  H20  molecule,  the  electronegativity  of  the  O  atomic  is  larger, 
so  that  the  H  atom  has  excess  force  to  form  a  hydrogen  bond  with 
highly  electronegative  atom  in  other  molecules,  such  as  F,  O  and  N. 
Fig.  4(a)  schematically  shows  the  hydrogen  bond  between 
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H 
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IV  A 
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Li 

Be 

B 

c 

N 
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F 
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Cl 
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3.0 

K 

Ca 

Ga 

Ge 

As 

Se 

Br 

0.8 

1.0 

1.6 

1.8 

2.0 
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2.8 

Rb 

Sr 

In 

Sn 

Sb 

Te 

1 
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1.0 
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1.8 
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2.5 
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Ti 

Pb 

Bi 

Po 

At 

0.7 

0.9 

1.8 

1.9 

1.9 

-- 

-- 

Fig.  3.  Some  atoms  electronegativity  values  (Pauling  scale). 


hydrophilic  IL  and  H20.  As  shown  in  [DMIM]C1  aqueous  solution, 
the  H  atom  in  H20  molecule  and  the  Cl  atom  in  IL  anion  of  form 
hydrogen  bond  0— H  ■  ■  •  Cl.  Fig.  4(b)  shows  the  hydrogen  bond¬ 
ing  interaction  between  R32  and  [BMlM]Tf2N  (l-butyl-3-methyli- 
midazolium  bis(trifluoromethylsulfonyl)imide).  The  imidazole  ring 
of  IL  has  three  hydrogen  bond  donors  (HI,  H2,  H3),  but  only  two 
proton  acceptors  (two  F  atoms)  in  the  R32  molecule,  so  two 
hydrogen  bonds  are  formed.  The  anion  [Tf2N]~  has  10  proton 
acceptors  (01,  02,  03,  04,  FI,  F2,  F3,  F4,  F5,  and  F6),  but  the  R32 
molecule  only  has  two  hydrogen  donors  (two  H  atoms),  so  two 
hydrogen  bonds  are  also  formed. 

Generally,  the  electronegativity  of  the  C  atomic  is  small,  and 
C— H  is  generally  difficult  to  form  hydrogen  bond.  Flowever,  the 
C  atom  between  two  N  atoms  in  the  imidazole  ring  of  cation  will 
form  a  hydrogen  bond  because  the  electronegativity  increases  due 
to  the  highly  electronegative  neighboring  two  N  atoms.  Therefore, 
in  IL  aqueous  solution,  the  0  atom  in  a  H20  molecule  can  easily 
form  a  hydrogen  bond  with  the  H  atom  on  the  C  atom  between 
two  N  atoms  in  the  imidazole  ring,  while  the  interaction  with  the 
H  atoms  on  the  other  two  C  atoms  in  the  imidazole  ring  is  weak. 
When  the  length  of  alkyl  side  chain  of  cation  in  IL  decreases,  the 
volume  of  cation  decreases,  the  surface  charge  density  increases, 
and  the  electrostatic  interaction  between  cation  and  H20  molecule 
enhances.  Correspondingly,  the  ability  of  IL  binding  together  with 
H20  molecule  enhances,  showing  that  a  shorter  cationic  alkyl  side 
chain  results  in  better  affinity  of  IL  with  H20.  Klimavicius  et  al.  [91  ] 
studied  hydrogen  bond  formed  in  D20  solutions  of  l-decyl-3- 
methylimidazolium  bromide  ([C10MIM]Br)  and  l-decyl-3- 
methylimidazolium  chloride  ([C10MIM]C1)  with  1H,  13C  NMR  and 
Raman  spectroscopy  monitoring  of  proton/deuteron  process.  They 
found  that  anions  and  aggregation  effects  play  crucial  role  in  IL 
aqueous  solutions,  anions  interact  with  cations  via  H-bond  and 
bind  water  molecules  in  their  solvation  shells. 

In  the  study  of  D'Angelo  et  al.  [92],  the  interaction  of  H20  and 
Br~  in  [BMIM]Br  was  investigated  combining  extended  X-ray 
absorption  fine  structure  spectroscopy  and  molecular  dynamics 
simulations.  A  picture  of  structural  properties  of  [BMIM]Br/H20 
was  gained,  in  the  surroundings  of  both  the  imidazolium  cation 
and  anion  Br~.  This  new  approach  has  been  found  to  be  particu¬ 
larly  well  suited  to  obtain  detailed  information  on  IL  solutions 
despite  their  complexity. 

For  the  working  pairs  adopting  water  as  refrigerant,  the 
interaction  between  water  and  absorbent  needs  to  be  enhanced, 
so  the  IL,  which  has  good  affinity  with  water,  is  studied.  The 
number  of  ILs  is  beyond  count,  because  the  structure  can  be 
designed.  The  different  configurations  of  the  anion  and  cation 
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Fig.  4.  Hydrogen  bonding  interaction  between  refrigerants  and  ILs.  (a)  H20/[DMIM]C1  and  (b)  R32/[BMIM]Tf2N. 


Table  3 

Relationship  of  1L  anion  and  water  solubility. 


Hydrophobic 

- 

Hydrophilic 

bf4- 

Br-,  Cl  ,  r 

[(CF3S02)2N] 

CF3S03- 

CH3COO-f  (CH3)2PC>4 

i(C2F5S02)2N]- 

CH30S03",  C2H5OSO3- 

make  ILs  show  different  hydrophilicities.  The  interaction  between 
IL  and  water  is  influenced  by  the  structure,  including  the  organic 
cationic  species,  cationic  substituent  structure,  and  anionic  spe¬ 
cies.  Studies  show  that  for  imidazole  ILs,  increasing  of  the  alkyl 
side  chain  of  the  cation  leads  to  reduce  the  hydrophilicities  of  ILs, 
such  as  [DMIM]X>[EMIM]X>[BMIM]X>[HMIM]X  (X  means  a 
certain  kind  of  anion). 

The  anion  characteristics  greatly  affect  the  solubility  of  IL  in 
water.  Table  3  lists  the  relationship  between  the  IL  anion  and 
water  solubility.  Generally,  for  the  solubility  of  the  refrigerant  in  IL, 
the  anion  shows  greater  effect  than  the  cation  [93],  because  the 
imidazolium  cation  only  depends  on  the  length  of  the  alkyl  side 
chain  or  structure  to  show  the  impact,  while  the  anion  may  have 
greater  variety  and  structural  changes. 

Similar  to  H20/IL  systems,  the  hydrogen  bond  has  a  very 
important  function  in  HFC/IL  systems.  The  hydrogen  bond  heavily 
influences  the  physical  and  chemical  properties  of  the  IL  working 
pairs.  Forming  hydrogen  bond  between  solute  and  solvent  will 
increase  the  solubility,  so  IL  having  good  affinity  with  HFC  is 
sought  for  HFC/IL  working  pairs.  As  in  hydrophilicities,  the  affinity 
of  IL  with  fluoride  is  also  influenced  by  the  cation  species, 
structure  of  cation  substitution,  and  anion  species.  Available 
experimental  results  show  that  for  imidazolium  ILs,  the  longer 
the  cationic  alkyl  side  chain  is,  the  better  the  affinity  with  HFC  will 
be,  as  shown  in  the  following  rule:  [EM1M]X<  [BMIM]X<  [HMIM] 
X<[OMIM]X  ([OMIM]  means  l-octyl-3-methylimidazolium), 
which  is  opposite  to  the  hydrophilic  of  IL  [75], 

Comparing  VLE  data  of  several  HFC/IL  systems,  the  effect  of 
anion  on  the  solubility  of  HFC  in  IL  was  as  follows:  [CnMIM]Tf2N  > 
[CnMIM]PF6  >  [CnMIM]TfO  ®  [CnMIM]BF4  [75],  The  solubility 


increases  with  increasing  number  of  fluorine  atoms  in  the  anion. 
The  hydrogen  bond  is  only  one  aspect  of  the  contribution  of  the  IL 
affinity  with  fluoride,  far  from  being  able  to  explain  the  pros  and 
cons  on  the  overall  performance  of  working  pairs.  Thus,  further 
studies  are  needed. 

2.3.  Infinity  dilution  activity  coefficient  criterion 

The  activity  coefficient  of  the  refrigerant  in  working  pairs 
relates  to  the  composition  of  the  system.  However,  the  infinite 
dilution  activity  coefficient  of  the  refrigerant  (yT),  calculated  on 
the  basis  of  Henry’s  law,  is  a  physical  quantity  that  is  unrelated  to 
the  composition  of  the  system  and  reflects  the  regular  interaction 
between  refrigerant  and  IL.  At  given  temperature  and  pressure, 
Henry’s  constant  in  working  pairs  is  given  by 

kt  =  Jbn  (p/xi)  =  yfPi(T)  (5) 

where  /q  is  Henry's  constant.  At  low  vapor  pressure  of  the 
refrigerant,  the  solubility  of  the  refrigerant  has  a  linear  relation¬ 
ship  with  vapor  pressure,  and  yf3  and  /q  can  be  obtained  when  the 
concentration  tends  to  zero,  i.e.,  the  slope  of  the  line  when 
pressure  tends  to  zero.  In  fact,  yf°  and  /q  are  direct  reflections  of 
the  solubility  of  the  refrigerant  at  low  vapor  pressure.  Smaller 
values  of  yf3  and  lq  mean  higher  solubility  of  the  refrigerant. 
Therefore,  yf3  and  Jq  can  be  used  to  compare  the  degree  of 
miscibility  of  the  refrigerant  and  IL. 

The  absorption  potential  was  proposed  as  another  evaluation 
criterion  for  absorption  cycle  working  pairs  [83], 

i//,  =  \/yf  (6) 

where  yq  is  the  absorption  potential.  At  given  temperature  and 
pressure,  higher  values  of  i/q  result  in  stronger  absorbing  ability  of 
the  absorbent. 

In  Table  4,  yf3  decreases  with  increasing  length  of  cation  alkyl 
side  chain.  In  Table  5,  except  for  l-butyl-3-methylimidazolium 
trifluoromethanesulfonate  ([BMIM]TfO),  yf3  decreases  with 
increasing  number  of  fluorine  atoms  of  anion.  The  yf3  of  binary 
systems  composed  of  nine  HFCs,  such  as  R32,  R134,  with  IL  [BMIM] 
PF6  are  listed  in  Table  6.  The  HFC  containing  more  carbon  atoms 
and  fluorine  atoms  shows  better  miscibility  with  [BMIM]PF6.  In 
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Table  4 

Activity  coefficients  at  infinite  dilution  of  R32  in  [PF6]-based  ILs  at  298.2  K. 


Item 

[OMIM]PF6 

|HMIM1PF6 

|BMIM]PF6 

[EMIM[PF6 

Absorbent 

structure 

\ _ /  F 

F^I^F 

F 

\  /\  /\  / 
N  N  +  v  v 

/PC 

f'pf 

F 

X  X]XXX/X 

\=/  F^  1  F 

/PC' 

F  I  F 

F 

\  ^  /X 

N  N  +  N  F 

\ _ /  FXLF 

F/PF 

F 

00 

n 

0.4525 

0.5559 

0.7328 

1.1028 

n 

2.2099 

1.7989 

1.3646 

0.9068 

Pls  /MPa 

1.6919 

1.6919 

1.6919 

1.6919 

k\ 

0.7656 

0.9405 

1.2398 

1.8658 

Table  5 

Activity  coefficients  at  infinite  dilution  of  R32  in  the  [BMlM]-based  ILs  at  298.2  K. 

Item 

[BMIM]Tf2N 

[BMIM]  TfO 

[BMIM[PF6 

[BMIM[BF4 

Absorbent 

AA  °  n  -? 

—  sx 

fA  O  o  Af 

F  F 

N  N  +  X'  x 

W  rF\ 

„  V*W\ 

°  W  a’ a 

W  r  -r 

structure 

F 

n  ~  F^t^F 

O  F 

1 

F 

co 

n 

0.5026 

0.7112 

0.7328 

0.8794 

1.9897 

1.4061 

1.3646 

1.1371 

p\  /MPa 

1.6919 

1.6919 

1.6919 

1.6919 

ki 

0.8503 

1.2033 

1.2398 

1.4879 

the  conductivity  will 

increase  by  10-20 

Table  6 

Activity  coefficients  at  infinite  dilution  of  various  HFCs  in  the  IL  [BMIM]PF6  at 
298.2  K. 


Items 

R32 

R134 

R41 

R23 

R152a 

Refrigerant  structure 

ch2f2 

chf2chf2 

ch3f 

chf3 

ch3chf2 

yT 

0.7328 

0.8289 

0.8677 

1.2847 

1.3957 

m 

1.3646 

1.2064 

1.1525 

0.7784 

0.7165 

ps,  /MPa 

1.6919 

0.5268 

3.8383 

4.7044 

0.5973 

Ri 

1.2398 

0.4367 

3.3305 

6.0437 

0.8337 

Items 

R161 

R134a 

R125 

R143a 

Refrigerant  structure 

ch3ch2f 

cf3ch2f 

chf2cf3 

ch3cf3 

yT 

1.4969 

1.8615 

3.1726 

3.6674 

m 

0.6680 

0.5372 

0.3152 

0.2727 

ps,  /MPa 

0.9230 

0.6664 

1.3797 

1.2632 

Ri 

1.3817 

1.2404 

4.3772 

4.6327 

association  between  the  anion  and  the  cation,  and  forming 
electroneutral  ionic  clusters  [95].  The  VLE  and  LLE  data  of  IL 
systems  can  be  well  described  by  the  CE  model  applied  for 
traditional  non-electrolyte  solutions,  such  as  the  NRTL  activity 
coefficient  model  and  PR  equation  of  state  (EOS),  showing  that  ILs 
can  be  reasonably  regarded  as  molecular  components  and  that  the 
electrostatic  interaction  can  be  ignored  [54,96],  Using  electrospray 
ionization  mass  spectrometry,  a  large  amount  of  ionic  clusters  in 
different  degrees  of  association  or  supramolecular  polymerization 
has  been  detected  in  ILs,  showing  that  ILs  are  similar  to  neutral 
molecules  or  weak  electrolytes  that  can  be  approximately  treated 
as  neutral  components. 

In  the  UNIFAC  model,  the  activity  coefficient  is  calculated  as  a 
sum  of  two  terms  [97], 


addition,  comparing  R134  with  R134a,  Y\  of  R134a  is  much  larger 
than  that  of  R134  because  of  isomerization. 

2.4.  VIE  prediction  without  experimental  data  by  the  UNIFAC  model 

The  activity  coefficient  of  the  refrigerant  species  in  IL  working 
pairs  is  the  basic  data  needed  for  selecting  suitable  IL  absorbents. 
However,  the  obtained  experimental  data  are  sometimes  insuffi¬ 
cient.  The  UNIFAC  model  provides  a  calculation  method  to  corre¬ 
late  and  predict  the  activity  coefficient  without  experimental  data, 
which  was  proposed  by  Fredenslund  based  on  the  group- 
contribution  method  in  1975  [94].  The  UNIFAC  model  is  specifi¬ 
cally  suitable  for  mixtures  of  molecules,  but  not  for  mixtures 
containing  ionic  components,  because  the  residual  term  of  the 
model  only  considers  short-range  interactions.  To  extend  the 
application  range  of  the  UNIFAC  model,  the  IL  can  be  considered 
as  a  neutral  molecule  or  a  weak  electrolyte  solution.  The  reason- 
ability  of  this  perspective  can  be  explained  from  the  following 
aspects  [81],  The  conductivity  of  IL  is  usually  very  low,  at 
approximately  10~3  mS/cm.  After  dilution  in  the  polar  solvent 


In  yt  =  In  j'f  +  In  yf  (7) 

where  yt  is  the  activity  coefficient  of  component  i,  yf  is  the 
combinatorial  term,  and  yf  is  the  residual  term.  The  combinatorial 
term  mainly  considers  the  difference  of  molecular  size  and  shape. 
The  residual  term  mainly  considers  the  interaction  force  between 
molecules,  and  can  be  obtained  from  the  properties  of  pure 
substances  and  mixtures. 

To  use  the  UNIFAC  model,  ILs  need  to  be  segmented  into 
groups.  ILs  can  be  decomposed  through  the  method  presented 
by  Kim  [98]  and  Lei  [99],  As  shown  in  Fig.  5,  the  IL  [BMIMjCl  is 
divided  into  one  CH3  group,  three  CH2  groups,  and  one  [M1M]C1 
group.  Other  ILs  can  also  be  divided  into  groups  in  the  same  way. 

Based  on  the  study  conducted  by  Lei  et  al.  [99],  Dong  et  al.  [71] 
studied  a  series  of  H20/IL  binary  systems,  which  mainly  included 
three  kinds  of  short  carbon  chain  imidazolium  cations,  namely, 
[DM1M],  [EMIM]  and  [BMIM],  as  well  as  seven  kinds  of  anions, 
namely,  Cl,  Br,  BF4,  DMP,  DBP  (dibutylphosphate),  EtS04,  and  TfO. 
Parameters  for  14  group  pairs  were  then  added  to  the  UNIFAC 
parameter  table.  The  proposed  new  group  interaction  parameters 
are  listed  in  Table  7,  which  can  be  used  to  predict  VLE  and  excess 
properties  of  systems  containing  H20  and  imidazolium  ILs  at 
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Fig.  5.  Group  segmentation  of  ionic  liquid  [BM1M]C1. 


various  compositions  and  temperatures.  The  new  parameters  are 
suitable  for  ILs  based  on  seven  kinds  of  anions,  namely,  Cl,  Br,  BF4, 
DMP,  DBP,  EtS04,  and  TfO.  Similarly,  works  have  been  carried  out 
for  HFC/1L  binary  systems  [75]  and  NH3/IL  binary  systems  [79],  The 
new  group  interaction  parameters  are  also  listed  in  Table  7,  which 
can  be  used  to  extend  the  evaluating  and  predicting  range  of  HFC/ 
IL  and  NH3/IL  systems  for  developing  alternative  working  pairs. 
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2.5.  Excess  property  criterion 

Closely  associated  with  yu  there  is  the  summation  relationship 
of  CE  [84], 

GE=  X  Xi[d(nGE)/dni]  T  =RT  £  xi ln  Yi  (8) 

i  =  1  i  =  1 

For  the  refrigerant  and  IL  binary  systems,  Eq.  (8)  can  be  written 
as 

GE  =  KT(x,  ln  Y\  +x2  In  y2)  (9) 

where  values  of  ln  jq  and  In  y2  can  be  determined  by  Eqs.  (l)-(3) 
using  experimental  data  or  with  the  help  of  the  UNIFAC  model. 

According  to  an  analysis  of  the  Flory-Huggins  equation  [101], 
when  the  attraction  of  heterogeneous  molecules  is  stronger  than 
that  of  homogeneous  molecules,  the  value  of  GE  is  less  than  zero. 
In  other  words,  a  negative  value  of  GE  means  the  affinity  between 
the  species  is  relatively  stronger. 

Thus,  for  a  binary  system  at  given  temperature  and  pressure,  a 
negative  GE  is  beneficial  for  solute  absorption,  and  the  smaller  the 
value  of  Ge  is,  the  better  the  absorption  effect  will  be.  At  given 
temperature  and  pressure,  GE  is  a  continuous  function  of  system 
composition,  which  may  exist  more  than  one  extreme  point.  The 
extreme  point  of  GE  presents  the  characteristic  of  thermodynamic 
properties  of  the  system.  Therefore,  the  extreme  point  of  GE  is 
presented  as  another  selection  criterion  for  evaluating  suitable 
absorbents  [102,103],  Thus,  at  given  temperature  and  pressure,  a 
smaller  value  of  the  extreme  point  of  GE  for  the  solute-absorbent 
binary  system  indicates  a  stronger  absorbing  ability  of  the 
absorbent. 

To  evaluate  the  solute-absorbent  binary  system,  the  criterion 
based  on  the  analysis  of  the  extreme  point  of  GE  is  called  “the 
criterion  of  minimum  extreme  value  of  the  GE-x  relationship”. 
The  criterion  based  on  the  analysis  of  the  absorption  potential  of 
solute  is  called  “the  criterion  of  maximum  value  of  the  absorption 
potential”.  The  former  is  a  macroscopic  trend  analysis  on  the  basis 
of  relatively  more  data  points,  while  the  latter  is  a  local  state  view 
that  requires  only  a  small  amount  of  data.  In  the  case  pressed  for 
available  VLE  data,  the  latter  may  be  more  convenient.  However, 
the  combination  of  two  methods  can  be  used  to  obtain  more 
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[MIM]Tf2N  145.80“  47.15  -100.69  0.0 

NH3  914.100  2.849  320.400  -  45.550  0.0  -  216.100 

[MIMJAc  1244.000  -388.000  0.0 
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comprehensive  evaluation  conclusions.  In  addition,  for  the  criter¬ 
ion  of  minimum  extreme  value  of  the  GE-x  relationship,  the 
absolute  value  of  minimum  extreme  value  of  the  CE  is  commonly 
marked  as  G|^ax. 

As  shown  in  Table  8  and  Fig.  6,  for  refrigerant  R32  absorption 
systems  with  absorbents  [BMIM]PF6,  [BMIM]BF4,  [BMIM]TfO,  and 
[EMIM]Tf2N,  the  Gj^ax  of  R32/[EMIM]Tf2N  binary  system  is  at 
minimum,  showing  that  [EMIM]Tf2N  has  better  ability  to  absorb 
the  refrigerant  R32. 

Comparing  data  from  the  UNIFAC  model  with  the  experimental 
data  listed  in  Table  8,  the  predicted  trends  of  G^ax  are  basically  the 
same  as  the  experimental  trends  of  the  y/x,  obtained  by  literature 
data.  The  assessment  results  obtained  by  the  G^ax  of  the  system 
agree  essentially  with  the  experimental  data  analysis.  For  R32 
species,  the  phase  equilibrium  ratio  values  of  the  R32/[EMIM]Tf2N 
system  is  relatively  smaller,  demonstrating  [EMIM]Tf2N  has  better 
absorbing  ability  compared  with  other  absorbents. 


3.  Researches  of  H20  and  IL  systems 

3.1.  Vapor  pressure 

As  presented  above,  excellent  working  pairs  are  generally  non¬ 
ideal  mixtures.  The  strong  affinity  of  absorbent  with  refrigerant 
often  leads  to  the  lower  saturated  vapor  pressure  of  solution  than 
that  of  ideal  solution.  In  other  words,  VLE  behavior  exhibits  highly 
negative  deviation  from  Raoult's  law,  which  is  used  to  describe  the 
VLE  behavior  of  ideal  mixture.  Therefore,  the  decrease  in  saturated 
vapor  pressure  of  solution  becomes  a  significant  indicator  for 
assessing  the  affinity  between  refrigerant  and  absorbent. 

In  recent  years,  many  researchers  have  investigated  the  VLE  of 
IL  working  pairs  systems.  Table  9  summarizes  the  VLE  measure¬ 
ment  of  H20/IL  systems  in  literatures. 

For  the  VLE  measurement  of  H20/IL  systems,  the  most  widely 
used  method  is  the  boiling  point  method.  The  principle  is  to 
precisely  prepare  a  mixture  with  known  composition,  observe 
phase  behavior  in  the  equilibrium  cell,  determining  the  properties 
in  equilibrium  state,  e.g.,  pressure  and  temperature. 

In  2004,  Kim  et  al.  [43]  measured  the  vapor  pressure  of  the 
systems  H20/[BMIM]Br,  H20/[BMIM]BF4,  and  H20/[HydeMIM]BF4 
( 1  -(2-hydroxyethyl)-3-methylimidazolium  tetrafluoroborate).  They 
proposed  that  these  systems  could  be  used  in  absorption  heat  pumps. 

Experimental  data  showed  that  the  vapor  pressure  of  water 
decreased  with  IL  addition,  but  the  decreasing  degree  differed 
according  to  IL  type  and  concentration.  For  imidazolium  based  ILs, 
the  interaction  with  H20  was  influenced  by  the  type  of  anion  and 
the  length  of  cationic  alkyl  side  chain.  Generally,  a  shorter  cationic 
alkyl  side  chain  leads  to  stronger  hydrophilicity.  Moreover,  the 
influence  of  the  anion  was  much  more  significant  than  that  of  the 
cationic  side  chain. 

In  recent  years,  researchers  have  studied  the  thermophysical 
properties  of  H20/phosphoric  IL  working  pairs,  such  as  [DMIM] 
DMP,  [EMIM]DMP,  l-ethyl-3-methylimidazolium  diethylphosphate 
([EMIM]DEP),  l-ethyl-3-ethyl-imidazolium  diethylphosphate  QEEIM] 


DEP),  and  l-butyl-3-methylimidazolium  dibutylphosphate  ([BMIM] 
DBP).  All  these  ILs  have  strong  affinity  with  water,  especially  [DMIM] 
DMP.  Dong  [72],  He  [109],  Kato  [112],  Wang  [117],  Zhao  [118]  et  al. 
measured  the  vapor  pressure  of  the  H20/[DMIM[DMP  system  at 
different  temperatures,  pressures,  and  compositions.  The  data  pro¬ 
vide  the  foundation  for  calculating  the  absorption  cycle,  and  con¬ 
firming  the  potential  of  the  H20/[DMIM[DMP  system  to  be  used  as 
working  pairs. 

Based  on  the  potential  use  of  hydrophilic  IL  as  absorbent 
species  combined  with  H20  to  constitute  alternative  working 
pairs,  Li  et  al.  [76]  chose  [DMIMJC1  and  l-methyl-3-methyl- 
imidazolium  tetrafluoroborate  ([DMIM]BF4)  as  additives,  added 
them  to  traditional  working  pairs  H20/LiBr  and  H20/LiCl,  and 
explored  the  modification  effect  of  IL  on  traditional  working  pairs. 
They  determined  the  saturated  vapor  pressures  of  ternary  systems 
at  different  temperatures  and  concentration  ranges  using  the 
boiling  point  method. 

The  correlation  method  of  VLE  data  generally  includes  two 
categories.  First  is  methods  only  using  equations  of  state,  such  as 
PR,  SRK,  BWR,  MH.  Second  is  methods  simultaneously  using 
equations  of  state  and  the  activity  coefficient  model,  such  as 
models  of  Margules,  Wilson,  NRTL,  UNIQUAC,  UNIFAC. 

For  H20/IL  systems,  Doker  and  Gmehling  [107]  analyzed  VLE 
data  correlation  of  H20/[EMIM[Tf2N  and  H20/[BMIM]Tf2N  systems 
using  three  kinds  of  activity  coefficient  models  including  Wilson, 
NRTL,  and  UNIQUAC.  Results  showed  that  NRTL  and  UNIQUAC 
models,  commonly  used  to  correlate  non-electrolyte  solutions,  are 
also  suitable  for  IL  working  pairs  with  good  accuracy.  Carvalho 
et  al.  [120]  determined  VLE  data  of  H20  combined  respectively 
with  l-ethyl-3-methylimidazolium  chloride  ([EMIM]C1),  [BMIM]C1, 
l-hexyl-3-methylimidazolium  chloride  ([HMIMjCl),  and  Choline 
Chloride  ([Nni(20H)]Cl),  the  NRTL  model  was  also  used  to  correlate 
with  data  and  a  good  fit  was  obtained. 

In  addition,  for  binary  and  ternary  systems  containing  H20  and 
IL,  Wang  et  al.  [117]  and  Li  et  al.  [76]  adopted  the  Antoine  equation 
to  correlate  VLE  data,  the  equation  is  expressed  as 

4 

log(p/kPa)  =  2  [Ai  +  lOOOBi/fT/K-dS.lSifw1  (10) 

i  =  0 


Fig.  6.  Ce-xr32  for  R32  systems  at  298.20  K.  ■,  [BMIM]BF4;  □,  [BMIMJTfO;  *, 
[BM1M]PF6;  A,  [EMIM]Tf2N. 


Table  S 

The  for  refrigerant  R32  and  IL  binary  systems  at  298.20  K. 


Refrigerant 

IL 

Experimental  data 

Cal.  by  UNIFAC 

r/K 

PI  MPa 

yjxi 

Ref. 

Gmax/J'  mol  1 

R32 

[EMIM]Tf2N 

C8H,,F6N304S2 

298.15 

0.399 

3.077 

[53] 

-408.47 

[BMIM]PF6 

c8h15f6n2p 

298.20 

0.400 

3.367 

154] 

-250.77 

[BMIM]BF4 

c8h15bf4n2 

298.20 

0.400 

3.676 

[54] 

-93.51 

[BMIMjTfO 

c9h15f3n2o3s 

298.14 

0.400 

3.788 

[73] 

-207.90 
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Table  9 


Summary  of  vapor  pressure  measurement  for  H20/IL  systems. 

IL 

Measuring  method 

Ref. 

Range  of  data 

T/K 

p/kPa 

[BMIMJCI 

Dynamic  recirculating  still  apparatus 

[104] 

373.15-441.11 

101.3 

[EMlM]EtS04 

Dynamic  recirculating  still  apparatus 

[105] 

373.15-410.58 

101.3 

[BMIM]MS04 

Dynamic  recirculating  still  apparatus 

[106] 

373.15-374.12 

101.3 

[EMlM]Tf2N 

Computer-driven  static  apparatus 

[107] 

353.15 

3.193-47.422 

[BMIM]Tf2N 

353.15 

5.317-47.156 

[DMlMjDMP 

Boilling  point  method 

[72] 

302.57-441.87 

2.14-101.01 

(EMIMJAc 

Static  apparatus 

[108] 

283.1-403.17 

0.215-104.8 

[HMIMJC1 

283.15-423.16 

0.298-262.8 

[DMlMjDMP 

Boilling  point  method 

[109] 

329.35-429.75 

6.62-39.15 

[EE1MJDEP 

Quasi-static  ebulliometer  method 

[110] 

317.63-371.31 

9.313-57.075 

(EMIM]BF4 

Gravimetric  method 

[111] 

373.15-403.15 

2.5-80.2 

[BMIMJBr 

Boilling  point  method 

[43] 

304.8-457.4 

3.5-103.2 

(BMIM]BF4 

311.4-475.2 

3.5-102.6 

[HydeMIM]BF4 

325.8-464.8 

3.3-100.8 

[DMlMjDMP 

Computer-driven  static  apparatus 

[112] 

353.15 

0.09-47.82 

[EMlM]Tf2N 

353.15 

3.19-47.39 

[BMIM]Tf2N 

353.15 

5.32-47.16 

(EMIMJDMP 

Bubble  point  method 

[113] 

320.35-411.15 

6.27-43.56 

(EMIMjTfO 

Dynamic  recirculating  still  apparatus 

[114] 

323.3-363.3 

5.9-69.3 

[EMIMJTFA 

328.3-368.2 

6.3-84.4 

[EMIM]EtS04 

Transpiration  method 

[115] 

302.9-322.9 

0.345-12.141 

[EMlMjDMP 

Quasi-static  ebulliometer  method 

[116] 

323.76-367.26 

12.609-58.804 

[DMlMjDMP 

Quasi-static  ebulliometer  method 

[117] 

328.936-379.411 

13.834-101.270 

[DMIMJCI 

Boiling  point  method 

[80] 

287.15-437.45 

1.10-100.98 

[DM1M]BF4 

Boiling  point  method 

[82] 

312.25-403.60 

4.69-103.23 

[DMlMjDMP 

Quasi-static  ebulliometer  method 

[118] 

312.88-373.01 

6.937-79.824 

[EMIMJDEP 

324.01-373.26 

11.872-83.919 

[BMIMJDBP 

317.98-368.07 

9.184-81.630 

[EMlM]EtS04 

Boiling  point  method 

[119] 

320.35-373.35 

4.21-24.75 

[EMIMJCl 

Boiling  point  method 

[120] 

355.58-435.03 

50-100 

[BMIMJC1 

354.92-430.10 

[HMIMJC1 

354.64-405.30 

[Nih(20H)]C1 

355.11-411.12 

[COC2MOR][FAP] 

Dynamic  method 

[121] 

290.1-351.7 

100 

[COC2MOR][Tf2N] 

291.0-349.6 

[C0C2PIP[[FAP1 

301.2-349.3 

[COC2PIP][Tf2N] 

313.2-353.9 

[COC2PYR][FAP] 

307.4-351.6 

[COC2PYR][Tf2N] 

291.8-357.5 

[BMP1P][N(CN)2] 

237.64-273.15 

[BMPYR)[N(CN)2] 

236.26-273.15 

[BMPy][N(CN)2] 

244.19-300.30 

[BMPyjjSCN] 

238.03-273.15 

where  T  is  the  temperature, 

and  w  is  the  mass  fraction  of  the 

Table  10 

absorbent.  A,  and  B,  are  equation  parameters  obtained  by  regres- 

Summary  of  heat  capacity  measurement  for  H20/IL  systems. 

sing  VLE  data  of  the  system. 

System 

Measuring  method 

Ref. 

Range  of  T/K 

3.2.  Heat  capacity 

H20/[DMIM[DMP 

Calvet  heat  flow  calorimeter 

[72] 

303.15-353.15 

H20/[EMIM[EtS04 

DSC 

[122] 

283.15-343.15 

Heat  capacity  is  an  important  factor  affecting  the  coefficient  of 
performance  for  absorption  refrigeration  cycles.  A  smaller  heat 

H20/[EMIM]TfO 

H20/[EMIM[TFA 

H20/[EMIM[EtS04 

DSC 

[123] 

283.15-343.15 

283.15- 343.15 

293.15- 318.15 

capacity  benefits  heat  transfer,  reduces  energy  consumption,  and 

H20/[BM1M]CH3S04 

293.15-318.15 

improves  the  cycle  performance.  Measuring  the  constant  pressure 

H20/[EMIM]TfO 

293.15-318.15 

heat  capacity  enables  the  calculation  of  enthalpy,  entropy,  Gibbs 

H20/[BMlM)TfO 

293.15-318.15 

free  energy,  and  other  thermodynamic  functions.  To  discuss  the 

H20/[DMIM]DMP 

H20/[EMIM]EtS04 

Adiabatic  solution  calorimeter 

DSC 

[109] 

[124] 

298.15-323.15 

303.2-353.2 

relationship  between  the  thermodynamic  model  of  new  working 

H20/iEMIM]TfO 

303.2-353.2 

pairs  and  circulation  mechanism,  the  heat  capacity  of  new  work- 

H20/[BM1M]CH3S04 

DSC 

[125] 

303.2-353.2 

ing  pairs  for  absorption  cycle  must  be  studied. 

H20/[BMIM]TfO 

303.2-353.2 

Previous  studies  on  the 

heat  caDacitv  determination  of  IL 

H20/[BM1M)BF4 

DSC 

[126] 

278.15-333.15 

working  pairs  in  literatures  mainly  focused  on  water  or  alcohol 
with  IL  systems,  whereas  the  heat  capacity  of  HFC  or  NH3  with  IL 
systems  was  rarely  reported  because  of  restricted  test  conditions. 

H20/[EMIM[DMP 

H20/[EMIM[EtS04 

Adiabatic  solution  calorimeter 
Isoperibol  solution  calorimeter 

[113] 

[119] 

298.15- 323.15 

303.15- 323.15 

Table  10  summarizes  researches  on  the  heat  capacity  of  H20/1L 

systems  published  in  recent  years. 

Heat  capacity  data  in  literatures  are  mostly  correlated  with 

Heat  capacity  is  commonly  measured  by  differential  scanning 

temperature  by  a 

polynomial  equation.  A  polynomial  equation  is 

calorimetry  (DSC)  and  the  micro-calorimeter  method. 

then  used  to  correlate  the  relationship  of 

heat 

capacity  with 
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temperature  and  composition. 

Cp=  £(a,  +  bj  T)x‘  (11) 

i  =  0 

where  Cp  is  the  heat  capacity,  a,  and  b,  are  equation  parameters 
obtained  by  regressing  experimental  data. 

3.3.  Density 

As  an  important  part  of  the  fluid  PVT  data,  liquid  density  plays  a 
significant  role  in  thermodynamics  and  is  necessary  for  calculating 
other  physical  properties  and  studying  of  heat  and  mass  transfer.  In 
practice,  the  density  of  working  pairs  is  closely  related  to  temperature 
and  composition,  so  the  influence  of  the  density  change  of  a  solution 
cannot  be  ignored  when  analyzing  the  absorption  refrigeration  cycle. 

The  density  of  most  imidazolium  based  ILs  roughly  ranges 
within  1.1  - 1.6  g/cm3,  which  is  more  important  than  that  of  water 
in  practical  applications.  The  density  of  1L  mainly  depends  on  the 
type  of  anion  and  cation,  and  the  influence  of  anion  is  more 
obvious.  Generally,  a  larger  anion  size  leads  to  greater  IL  density, 
and  a  larger  volume  of  organic  cation  leads  to  smaller  IL  density. 

Previously  studies  on  the  density  determination  of  IL  working  pairs 
in  literatures  mainly  focused  on  water  or  alcohol  with  IL  systems, 
whereas  the  density  of  HFC  or  NH3  with  IL  systems  was  rarely 
reported  because  of  restricted  test  conditions.  Researches  on  the 
density  of  H20/IL  systems  in  recent  years  are  summarized  in  Table  11. 

Density  is  commonly  measured  by  the  pycnometer  [127,128], 
gravity  balance  (Westphal)  [109,119,129,130],  and  vibrating 
U-tube  densimeter  [74,  131  - 139].  The  vibrating  U-tube 


Table  11 

Summary  of  density  measurement  for  H20/IL  systems. 


System 

Measuring  method 

Ref. 

Range  of  T/K 

H20/[BMIMJN(CN)2 

Vibrating  tube  densimeter 

[131J 

278.15-363.15 

H20/[BMIM[C(CN)3 

278.15-363.15 

H20/[DMIM]C1 

Vibrating  tube  densimeter 

[74] 

293.15-318.15 

H20/[BMlM)Br 

Vibrating  tube  densimeter 

[332] 

298.15 

H20/[DMIM]Br 

298.15 

H20/[EMIM)EtS04 

Vibrating  tube  densimeter 

[123] 

293.15-318.15 

H20/[BMIM[CH3S04 

293.15-318.15 

H20/[EMIM]TfO 

293.15-318.15 

H20/[BMIM]TfO 

293.15-318.15 

H20/[BM1M]C1 

Vibrating  tube  densimeter 

[133] 

298.15-343.15 

H20/[HMIM]C1 

298.15-343.15 

H20/[0MIM]C1 

298.15-343.15 

H20/[EMIM)EtS04 

Vibrating  tube  densimeter 

[134] 

298.15-328.15 

H20/[BMIM]CH3S04 

Vibrating  tube  densimeter 

[135] 

298.15-328.15 

H20/[DMIM]DMP 

Gravity  balance 

[109] 

298.15-323.15 

H20/[BMIM]C1 

Density  bottle 

[136] 

298.15 

H20/[EMIM)Br 

298.15 

H20/[BMIM]Br 

298.15 

H20/[BMIM)BF4 

298.15 

H20/[EMIM]EtS04 

Westphal  balance 

[129] 

278.15-333.15 

H20/[EMIM]BF4 

Vibrating  tube  densimeter 

[137] 

298.15 

H20/[BMIM]BF4 

298.15 

H20/[HMIM]BF4 

298.15 

H20/[EMIM]EtS04 

Vibrating  tube  densimeter 

[138] 

278.15-348.15 

H20/[EMIM]TfO 

278.15-348.15 

H20/[EMIM]TFA 

278.15-348.15 

H20/[EMIM)BF4 

Vibrating  tube  densimeter 

[95] 

298.15 

H20/[BMIM)BF4 

298.15 

H20/[HMIM]BF4 

298.15 

H20/[BMIM]PF6 

298.15 

H20/[BMIM[DCA 

298.15 

H20/[EMIM]BF4 

Vibrating  tube  densimeter 

[139] 

298.15 

H20/[EMIM]EtS04 

298.15 

H20/[EMIM]EtS04 

Westphal  balance 

[130] 

278.2-338.2 

H20/[EMIMJBF4 

Pycnometer 

[127] 

293.15-323.15 

H20/[BMIM]BF4 

Pycnometer 

[128] 

303.15-353.15 

H20/[EMIM[EtS04 

Gravity  balance 

[119] 

303.15-323.15 

densimeter  method  is  one  of  the  most  common  and  highest 
accuracy  methods  to  date.  The  uncertainty  can  reach 
+ 1  x  10~6gcm~3,  and  the  measurement  error  does  not  exceed 
0.005%.  The  working  principle  is  to  change  the  inherent  oscillation 
frequency  of  the  vibrating  tube  with  the  mass  change  (i.e.,  the 
density  change)  of  the  liquid  flowing  through  the  vibrating  tube. 
After  determining  the  relation  between  the  oscillation  period  of 
the  vibrating  tube  and  liquid  density,  the  density  can  be  directly 
measured. 

Solution  density  can  be  correlated  with  composition  and 
temperature  by  various  expressions.  Given  that  the  density  of 
the  IL  solution  has  a  simple  relation  with  temperature  and 
composition,  a  polynomial  equation  is  selected  to  regress  the 
density  data  of  H20/IL  systems  as 

p=  £  A„wn  +  T  2  B„wn  +  T2  £  C„wn  (12) 

n  =  0  n  =  0  n  =  0 

where  p  is  the  density  of  the  binary  system,  w  is  the  mass  fraction 
of  water.  An,  Bn  and  Cn  are  polynomial  parameters  obtained  by 
fitting  experiment  data  through  the  least  square  method. 

3.4.  Other  properties 

For  a  detailed  investigation  of  the  suggested  working  pairs, 
more  properties  such  as  viscosity,  surface  tension,  and  thermal 
conductivity  are  required. 

The  viscosity  of  IL  working  pair  systems  was  rarely  reported  in 
literatures.  The  studied  systems  are  mainly  H20/[EMIM]BF4,  H20/ 
[BMIM]BF4,  H20/[HMIM]BF4i  H20/[BMIM]C1,  H20/[EMIM]Br  (1-ethyl- 
3-methylimidazolium  bromide),  H20/[BMIM]Br,  H20/[EMIM[EtS04 
(l-ethyl-3-methylimidazolium  ethylsulfate),  H20/[BMIM]MeS04  (1- 
butyl-3-methylimidazolium  methylsulfate),  and  so  on.  Results  show 
that  the  viscosity  of  H20/IL  systems  gradually  decreases  with  the 
increasing  water  content.  The  viscosity  of  IL  is  generally  tens  to 
hundreds  times  higher  than  that  of  water  mainly  because  of  the 
different  degrees  of  IL  aggregation,  due  to  the  strong  hydrogen  bond 
between  anion  and  cation.  When  water  is  added  to  IL,  the  hydrogen 
bond  between  the  anion  and  cation  weakens,  which  increases  the 
mobility  of  ions  and  decreases  the  viscosity  of  IL 

The  surface  tension  of  IL  is  between  that  of  water  and  common 
solvent,  and  decreases  with  the  increase  of  temperature.  However, 
the  surface  tensions  of  the  refrigerant/lL  systems  are  still  rarely 
reported.  The  systems  conducting  research  on  surface  tension 
include  H20/[EMIM]BF4,  H20/[BMIM]BF4,  H20/[HMIM]BF4,  H20/ 
[BMIMJC1,  H20/[BMIM]Br,  H20/[EMIM]Br,  and  so  on.  The  addition 
of  IL  can  reduce  the  surface  tension  of  the  refrigerant  quickly, 
which  is  beneficial  to  application  in  absorption  refrigeration  cycle. 

The  reports  on  the  thermal  conductivity  of  IL  working  pair  systems 
which  is  necessary  for  the  design  of  heat-transfer  equipment  are 
limited.  The  thermal  conductivities  for  the  binary  systems  H20 / 
[EMIM]EtS04,  H20/[DMIM]DMP,  CH3OH/[DMIM]DMP,  and  C2H5OH/ 
[EMIM]EtS04  were  reported.  It  is  expected  that  refrigerant/IL  systems 
could  be  perfect  heat-transfer  fluids. 

Moreover,  properties  such  as  toxicity,  environmental  impact,  and 
cost  have  to  be  evaluated  and  considered  for  further  investigation. 

4.  Researches  of  other  refrigerant  and  IL  working  pair  systems 

4.3.  NH3  and  IL  systems  (solubility  and  other  properties) 

Because  of  the  good  mutual  solubility  of  NH3  and  H20  and  the 
great  latent  heat  of  evaporation,  the  NH3/H20  working  pair  is 
widely  used  in  various  kinds  of  absorption  refrigeration  cycles. 
However,  the  system  has  some  drawbacks  such  as  high  working 
pressure,  toxicity,  and  difficulty  in  regenerating  NH3  from  H20. 
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Table  12 

Summary  of  solubility  measurement  of  NH3/IL  systems. 


Ionic  liquid 

Measuring  method 

Ref. 

Range  of  T/K 

Range  of  p/ MPa 

Range  of  x/mol% 

[BMIM]PF6 

283.4-355.8 

0.138-2.700 

23.9-86.2 

[HMIMJC1 

Static  method 

[231 

283.1-347.9 

0.044-2.490 

6.0-83.7 

[EMlM]Tf2N 

283.3-347.6 

0.114-2.860 

4.5-94.8 

|BMIM]BF4 

282.2-355.1 

0.091-2.570 

6.8-84.4 

|EMIM]Ac 

282.5-348.5 

0.321-2.891 

47.3-87.7 

[EMlM]EtS04 

Static  method 

[501 

282.7-372.3 

0.287-4.777 

42.4-87.5 

[EMlMjSCN 

283.2-372.8 

0.244-5.007 

34.0-87.6 

[DMEAJAc 

283.2-372.8 

0.136-4.249 

45.4-86.5 

[EM1M]BF4 

293.15-333.15 

0.11-0.63 

11.85-69.21 

|BMIM]BF4 

Static  method 

[140] 

293.15-333.15 

0.07-0.83 

6.08-75.31 

|HMIM]BF4 

293.15-333.15 

0.14-0.71 

12.80-75.43 

[0M1M]BF4 

293.15-333.15 

0.10-0.61 

13.21-80.81 

[DMIMJDMP 

Isothermal  synthetic 

[78] 

293.15-333.15 

0.053-0.665 

7.0-75.2 

[BMIM]Zn2Cl5 

Static  method 

[141] 

323.15-563.15 

0.067-1.965 

83.62-94.67 

Therefore,  a  great  number  of  new  absorbents  including  1L  have 
been  studied  as  alternatives. 

To  find  novel  working  pairs  containing  NH3  as  refrigerant  and  1L  as 
absorbent,  researchers  have  carried  out  related  works.  The  research 
progress  of  VLE  of  NH3/IL  systems  is  presented  in  Table  12.  The 
common  method  of  determining  the  VLE  data  of  NH3/IL  systems  is  as 
follows.  First  a  certain  mass  of  1L  is  added  to  the  static  equilibrium  cell 
that  is  then  filled  with  a  certain  mass  of  NH3.  Second  the  equilibrium 
cell  is  placed  in  a  thermoregulated  system  (oil  or  water  bath).  The 
temperature  and  pressure  are  recorded  when  the  NH3/IL  system 
reaches  equilibrium  at  constant  temperature.  Third  the  equilibrium 
pressures  at  different  temperatures  are  recorded  by  varying  the 
temperature.  Finally  an  appropriate  model  is  used  to  calculate  and 
obtain  p—T—x  data.  Using  these  experimental  data,  researchers  can 
then  verify  the  applicability  of  the  thermodynamic  model  to  NH3/IL 
working  pairs,  simulate  the  working  condition  of  the  absorption 
cycle,  and  explore  the  prospects  of  the  working  pairs  in  practical 
applications. 

As  shown  in  Table  12,  Yokozeki  et  al.  [23,50]  pioneered  works 
in  this  region.  They  found  that  solubility  of  NH3  in  ILs  such  as 
[HMIM]C1  and  [BM1M]PF6  at  298  K  is  almost  equivalent  to  that  in 
water,  i.e.,  negative  deviations  from  Raoult's  law  for  the  two 
systems  are  at  relatively  the  same  level.  Another  important  result 
is  that  NH3  does  not  react  with  ILs  under  the  experimental 
condition,  and  no  loss  of  solvent  occurs  even  in  the  acid  solvent 
[HMIM]C1.  However,  no  absorption  mechanism  between  NH3  and 
various  ILs  is  proposed  and  further  study  is  deserved  at  the 
molecular  level  by  spectroscopy  methods. 

In  another  work,  Yokozeki  [50]  calculated  residual  properties 
such  as  residual  enthalpy,  residual  entropy,  and  residual  Gibbs 
energy  of  NH3/IL  systems.  They  found  that  these  residual  proper¬ 
ties  are  all  negative.  The  negative  value  of  the  NH3/[DMEA]Ac 
system  was  the  maximum,  and  residual  properties  of  this  system 
were  close  to  those  of  the  NH3/H20  system,  indicates  that  a  strong 
interaction  like  hydrogen  bond  exists  between  NH3  and  ILs. 

In  2010,  Li  et  al.  [140]  determined  the  VLE  data  of  NH3  and  four 
ILs  [CnMIM]BF4  (n  =  2,  4,  6,  and  8)  and  studied  the  effect  of  alkyl 
side  chain  length  of  ILs  on  NH3  absorption.  Results  show  that  the 
solubility  of  NH3  in  studied  ILs  is  high  and  increases  with  the 
increase  of  alkyl  side  chain  length.  They  believed  that  increase  of 
carbon  chain  length  leads  to  the  decrease  of  ILs  density,  generating 
more  free  volume,  and  dissolving  more  NH3  molecules.  Using  the 
Krichevisky-Kasarnovsky  equation,  Li  et  al.  [140]  calculated  the 
Henry's  constant  and  partial  molar  volume  of  NH3i  they  estab¬ 
lished  a  thermodynamic  model  to  calculate  the  enthalpy,  entropy, 
Gibbs  energy,  and  heat  capacity  of  the  system.  Results  confirm  the 
non-spontaneity  of  the  absorption  process  and  the  disadvantage  of 
high  temperature  in  NH3  absorption. 


Chen  et  al.  [141]  measured  total  pressures  of  NH3  and  metal 
ion-containing  IL  [BMIM]Zn2Cl5  by  a  static  method,  and  correlated 
experimental  data  by  the  modified  UNIFAC  model.  The  solubilities 
of  NH3  in  [BMIM]Zn2Cl5  is  higher  than  that  in  the  other  ILs  tested 
by  Yokozeki  et  al.  [50],  but  lower  than  that  in  ZnCl2.  Author's  group 
[78]  chose  NH3  and  [DMIM]DMP  as  working  pairs  and  determined 
solubilities  by  the  isothermal  synthesis  method.  The  solubility 
decreases  with  increased  temperature  and  increases  with 
increased  pressure.  The  experimental  data  are  consistent  with 
the  calculated  results  from  the  NRTL  model,  with  a  relative 
deviation  within  5%. 

The  literature  survey  reveals  that  properties  of  NH3/IL  working 
pairs  have  been  gradually  examined  in  recent  years.  Given  the 
controllability  of  the  IL  structure  (i.e.,  various  ILs  can  be  synthe¬ 
sized  according  to  a  given  target  or  practical  needs),  the  gap  in  this 
applied  research  field  urgently  needs  to  be  filled.  Moreover, 
considering  the  high  solubility  and  low  energy  consumption  in 
the  separation  of  NH3/IL  systems,  it  is  necessary  to  study  the 
interaction  between  NH3  and  various  ILs,  for  the  development  of 
suitable  working  pairs  for  absorption  cycle. 

4.2.  HFCs  and  IL  systems  (solubility  and  other  properties) 

As  a  kind  of  refrigerant,  HFCs  are  also  concerned  in  developing 
novel  absorption  IL  working  pairs.  The  research  progress  of  the 
VLE  of  HFC/IL  systems  is  summarized  in  Table  13.  The  studied 
refrigerants  are  tetrafluoromethane  (R14),  trifluoromethane  (R23), 
R32,  fluoromethane  (R41),  R125,  R134,  R134a,  R152a,  and  fluor- 
oethane  (R161). 

Shiflett  and  Yokozeki  et  al.  [47,53,54,58,61,62]  determined  the 
solubility  and  diffusivity  of  HFCs  such  as  R32,  R134a,  R125,  R23, 
R152a,  R134,  and  R161  in  ILs  by  gravimetric  microbalance,  and 
used  the  EOS,  NRTL,  and  Antoine-type  equation  to  process  the 
experimental  data.  Ren  et  al.  [96]  determined  the  solubility  of 
R134a  in  [HMIM]-based  ILs,  regressed  the  experimental  data  using 
the  PR  EOS  with  van  der  Waals  two-parameter  mixing  rules,  and 
studied  the  effect  of  the  anion  and  cationic  carbon  chain  length  on 
the  solubility  and  molar  volume  of  R134a  in  ILs.  Author's  group 
[73]  determined  the  solubility  of  R32  and  R152a  in  l-ethyl-3- 
methylimidazolium  trifluoromethanesulfonate  ([EMIMjTfO)  and 
[BMIM]TfO  by  the  isothermal  synthesis  method,  results  well  fit 
the  NRTL  model. 

4.3.  HC  and  IL  systems  (solubility  and  other  properties) 

Many  researchers  have  studied  the  physical  and  chemical 
properties  of  HC  and  IL  binary  systems  to  discuss  their  potential 
use  as  working  pairs  of  absorption  cycle.  Lee  et  al.  [144]  measured 
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Table  13 

Summary  of  solubility  measurement  of  HFC/IL  systems. 


Refrigerant 

Absorbent 

Measuring  method 

Ref. 

Range  of  T  /K 

Range  of  p/ MPa 

R32 

[EMIMJTfO 

Isothermal  synthetic 

[73] 

273.14-298.15 

0.097-0.857 

[BMIMJTfO 

273.13-298.19 

0.104-0.902 

R152a 

[EMlM]TfO 

273.16-348.17 

0.040-0.848 

[BMIMJTfO 

273.22-348.15 

0.056-0.879 

R14 

[HMIM]Tf2N 

Synthetic  method 

[142] 

293.3-413.3 

1.202-9.582 

R134a 

[HMIM]Tf2N 

Static  method 

[96] 

398.15-348.15 

0.042-4.320 

[HMIM]PF6 

0.081-14.330 

[HMIM]PF6 

0.081-12.880 

R23 

[EM1M]PF6 

Synthetic  method 

[143] 

308.17-367.33 

1.602-51.64 

R32 

[BMIM]PF6 

Gravimetric  microbalance 

[54] 

283.2-348.2 

0.0097-0.9999 

[BMIM]BF4 

283.0-348.2 

0.0097-0.9999 

R125 

[BMIM]PF6 

283.1-348.3 

0.0099-0.9998 

R134a 

[BMIM]PF6 

283.0-348.2 

0.0097-0.3500 

R23 

[BMIM]PF6 

282.6-348.1 

0.0098-2.0002 

R152a 

[BMIM]PF6 

283.1-348.2 

0.0097-0.4505 

R32 

[DMPlM]TMeM 

Gravimetric  microbalance 

[53] 

283.15-348.15 

0.0094-1.0005 

[EMIMJBEI 

283.15-348.05 

0.0096-1.0005 

[DMPIMJBMel 

298.15 

0.0099-1.0011 

[EMIMJBMel 

283.15-348.05 

0.0096-1.0005 

[Pmpy]BMeI 

283.15-348.05 

0.0095-1.0004 

[bmpy]BMeI 

298.15 

0.0096-1.0000 

[BMIMJAc 

0.0099-1.0004 

[BMIMjSCN 

0.0095-0.9992 

[BMIM]MeS04 

0.0099-1.0006 

[EMIMJTFES 

0.0099-1.0016 

[BMIMJTFES 

0.0097-0.9989 

[HMIMJTFES 

0.0099-0.9980 

[DMIMJTFES 

0.0096-1.0010 

[BMIMJHFPS 

0.0095-1.0004 

[BMIMjFS 

0.0100-1.0005 

[BMIMjTPES 

0.0095-0.9994 

JBMIMjTTES 

0.0095-0.9992 

R134a 

[EM1MJBE1 

Gravimetric  microbalance 

[62] 

283.10-348.10 

0.0103-0.3505 

[BMIM]HFPS 

283.10-348.10 

0.0099-0.3506 

[BMIMjTPES 

283.05-348.10 

0.0102-0.3505 

[BMIMjTTES 

283.10-348.15 

0.0102-0.3505 

[6,6,6,14-P]TPES 

282.90-348.10 

0.0098-0.3504 

[4,4,4, 14-PJHFPS 

283.05-348.10 

0.0099-0.3504 

R41 

[BMIM]PF6 

Gravimetric  microbalance 

[58] 

283.09-348.18 

0.0097-1.9995 

R161 

[BMIM]PF6 

283.06-348.16 

0.0099-0.7005 

R134 

[BMIM]PF6 

283.11-348.16 

0.0101-0.3505 

R134 

[EMIM]Tf2N 

Gravimetric  microbalance 

[47] 

282.9-348.1 

0.0099-0.3505 

R125 

[EMIM]Tf2N 

Gravimetric  microbalance 

[61] 

283.1-348.2 

0.0100-0.9998 

R134a 

[EMIM]Tf2N 

282.7-348.1 

0.0100-0.3503 

solubility  data  of  various  HCs  in  [BMIM]Tf2N  by  the  saturation 
method  and  obtained  Henry’s  constant  through  the  experimental 
data.  Shiflett  et  al.  [47,60,61]  determined  the  solubility  data  of  a 
series  of  HCs  and  fluorinated  benzene  in  [EMlM]Tf2N  using  the 
gravimetric  microbalance  and  volumetric  method,  then  experi¬ 
ment  data  were  correlated  with  the  NRTL  model.  Liu  et  al. 
[145,146]  predicted  solubilities  of  small  HCs  in  about  700  ILs, 
and  reported  solubilities  of  methane,  ethane,  ethylene  and  pro¬ 
pane  in  trimethyloctylphosphonium  bis(2,4,4-trimethylpentyl) 
phosphinate  ([P8111]TMPP),  tetrabutylphosphonium  bis(2,4,4-tri- 
methylpentyl)phosphinate  ([P4444]TMPP)  and  [EMIM]Tf2N. 
Table  14  summarizes  studies  on  the  VLE  of  HC/IL  systems  pub¬ 
lished  in  recent  years. 

4.4.  Alcohol  and  1L  systems  (solubility  and  other  properties ) 

Kim  et  al.  [42]  determined  the  solubility  data  of  TFE  in 
[BMIM]Br  and  [BMIM]BF4  by  the  boiling  point  method  and 
correlated  the  data  with  the  Antoine-type  equation.  Shiflett 
et  al.  [56]  determined  the  VLE  data  of  the  n-butanol/[BMlMJPF6 
system  using  the  volumetric  method.  Verevkin  [148],  Zhao 
[118],  Wang  [116],  Jiang  [110],  Wang  [80,117],  Shen  [149]  and 
Carvalho  et  al.  [120]  determined  the  VLE  data  of  alcohols  and 
benzene  in  different  ILs  by  the  static  method,  correlated  the 


data  using  NRTL,  Antoine-type  and  other  equations,  and 
obtained  the  activity  coefficients  and  model  parameters  of 
different  solutes  in  ILs.  Gonzalez  et  al.  [150]  determined  vapor 
pressure  of  four  alcohols  in  l-hexyl-3-methylimidazolium  tri- 
fluoromethanesulfonate  ([HMIM]TfO)  using  the  vapor  pressure 
osmometry  technique  and  derived  other  thermodynamic  prop¬ 
erties.  Particularly,  Zhang  et  al.  [151]  proposed  a  new  alter¬ 
native  working  pairs  for  absorption  refrigeration  cycle,  which 
was  ternary  solutions  H20/CH3OH(or  C2H5OH)/[BMIM]DBP,  and 
measured  vapor  pressure  of  binary  solutions  H20(CH3OH, 
orC2H5OH)/[BMIM]DBP  with  an  inclined  boiling  apparatus  to 
obtain  binary  interaction  parameters  of  NRTL  activity  coeffi¬ 
cient  model,  then  measured  and  predicted  vapor  pressure  of 
proposed  solutions  with  the  same  apparatus  and  model.  Studies 
on  the  VLE,  density,  heat  capacity,  and  viscosity  of  alcohol  and 
IL  systems  published  in  recent  years  are  listed  in  Tables  15-18. 


5.  Assessment  of  absorption  cycle  adopting  IL  working  pairs 

5.3.  Simulation  for  single-effect  absorption  cooling  cycle 

To  determine  whether  a  working  pair  has  the  potential  for 
further  research  and  development,  it  should  be  placed  in  an 
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Table  14 

Summary  of  solubility  measurement  of  HC/IL  systems. 


Refrigerant 

Absorbent 

Measuring  method 

Ref. 

Range  of  T  /K 

Range  of  p  /kPa 

Propane 

[BMIM]Tf2N 

Saturation  method 

[144] 

279.98-339.97 

88.3-1062.2 

Propene 

999-1218.2 

Butane 

280.00-340.00 

334-304.1 

1 -Butane 

21.0-351.1 

Hexane 

[BMIMJTfO 

Static  apparatus 

[147] 

363.15 

15.79-165.7 

1 -Hexane 

22.54-180.9 

2,2,4-Trimethylpentane 

11.86-65.21 

1-Nonene 

2.599-16.34 

Decane 

2.000-6.202 

1 -Decane 

1.132-6.989 

1 -Hexane 

[OMlMJTfO 

16.19-198.8 

1-Octene 

1.672-30.24 

1-Nonene 

0.784-13.45 

Fluorinated  benzenes 

(EMlM]Tf2N 

Volumetric  method 

]60] 

282.5-373.0 

- 

Benzene 

[BMIM]Tf2N 

Static  method 

[148] 

298.15-313.15 

0-24.331 

Methane 

[P8111]TMPP 

Isochoric  saturation  method 

[145,146] 

299-323 

1090-3610 

Ethane 

540-1690 

Ethylene 

1080-3540 

Propane 

50-210 

Methane 

Ethane 

[P4444]TMPP 

313-353 

0-5000 

Ethylene 

Propane 

Methane 

Ethane 

Ethylene 

Propane 

[EMIM]Tf2N 

299-354 

0-4000 

Table  15 

Summary  of  vapor  pressure  measurement  of  alcohol/IL  systems. 

Refrigerant 

Absorbent 

Measuring  method 

Ref. 

Range  of  T  /K 

Range  of  p/kPa 

Ethanol 

[EMIMJTOS 

Ebulliometric  method 

[152] 

373.15 

0.0533-0.2250 

Methanol 

[EEIMJDEP 

Quasi-static  method 

[110] 

299.20-366.19 

14.336-65.452 

Ethanol 

[EEIMJDEP 

309.80-370.80 

9.139-56.716 

Methanol 

[EMIM]DMP 

Bubble  point  method 

[113] 

301.95-386.15 

0.01134-0.04630 

Ethanol 

309.05-389.15 

0.00747  -  0.05110 

Methanol 

[DMIMJC1 

Quasi-static  method 

[149] 

298.85-348.63 

16.65-102.57 

Ethanol 

[DMIMJCl 

312.19-361.04 

11.54-101.90 

1 -Butanol 

[BMIM]PF6 

Volumetric  method 

[56] 

285.7-332.0 

- 

Methanol 

[BMIM]Tf2N 

Static  method 

[148] 

298.15-313.15 

0-35.450 

Ethanol 

[BMIM]Tf2N 

0-17.928 

Propanol 

[BMIM]Tf2N 

0-6.986 

Methanol 

[EMIMJDMP 

Quasi-static  method 

[116] 

299.20-343.86 

15.969-60.419 

Ethanol 

[EMIM]DMP 

309.95-350.67 

14.165-60.073 

1 -Propanol 

[DMIMJDMP 

Quasi-static  method 

[117] 

330.478-401.916 

11.094-100.935 

2-Propanol 

[DMIMJDMP 

316.402-384.177 

11.149-102.082 

TFE 

[BMIMJBr 

Boiling-point  method 

[42] 

293.2-469.4 

0.0049-0.1010 

[BMIM]BF4 

315.2-469.4 

0.0068-0.1003 

[EMIM]BF4 

Boiling-point  method 

[80] 

283.15-460.85 

0.00306-0.10102 

Methanol 

[DMIMJDMP 

Quasi-static  method 

[118] 

285.96-343.71 

7.811-76.224 

[EMIMJDEP 

285.31-342.19 

7.943-79.563 

[BMIMJDBP 

286.82-335.10 

8.431-80.194 

Ethanol 

[DMIMJDMP 

298.62-359.57 

7.853-74.249 

[EMIMJDEP 

302.30-356.00 

9.068-81.091 

[BMIMJDBP 

298.42-355.12 

7.606-77.658 

Ethanol 

[EMIM]C1 

Boilling  point  method 

[120] 

334.80-392.42 

50-100 

[BMIMJC1 

335.03-405.89 

[HMIMJC1 

335.15-411.86 

[Nih(20H)]C1 

334.99-353.06 

h2o/ch3oh 

[BMIMJDBP 

Boiling-point  method 

[151] 

315.95-392.65 

11.34-72.02 

h2o/c2h5oh 

321.65-395.65 

11.76-69.02 

1 -Propanol 

[HMIMJTfO 

Vapor  pressure  osmometry  technique 

[150] 

323.15 

10.39-12.07 

2-Propanol 

20.12-23.52 

1 -Butanol 

3.89-4.46 

2-Butanol 

8.85-10.56 

1-Pentanol 

1.44-1.63 

absorption  cycle  to  assess  the  performance  by  simulation  or 
experiment.  Some  researchers  have  carried  out  such  works,  but 
the  main  method  is  still  limited  to  simulation  and  analysis  on 


the  basis  of  thermophysical  properties  of  the  IL  working 
pairs.  The  works  mainly  focus  on  single-effect  absorption 
cooling  cycle. 
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Table  16 

Summary  of  density  measurement  of  alcohol/IL  systems. 


Refrigerant 

Absorbent 

Measuring  method 

Ref. 

Range  of  T/K 

Ethanol 

Methanol 

[DMIMJDMP 

Gravity  balance 

[109] 

298.15-323.15 

Methanol 

[OMIMJC1 

Vibrating  tube  densimeter 

[153] 

298.15-328.15 

Methanol 

[EMIMJDMP 

[BMIMJDMP 

Vibrating  tube  densimeter 

[154] 

293.15-333.15 

Ethanol 

[EMIMJDMP 

[BMIMJDMP 

Ethanol 

[EMIMJHS 

[EMIMJOS 

[EMIMJES 

[EMIMJBS 

Viscodensimeter 

[155] 

288-318 

Ethanol 

[EMIM]Tf2N 

Density  meter 

[156] 

293.15-323.15 

Table  17 

Summary  of  heat  capacity  measurement  of  alcohol/IL  systems. 

Refrigerant 

Absorbent 

Measuring  method 

Ref. 

Range  of  T/K 

Ethanol 

Methanol 

[DMIMJDMP 

Adiabatic  solution  calorimeter 

[109] 

298.15-323.15 

Ethanol 

Methanol 

[EMIMJDMP 

Adiabatic  solution  calorimeter 

[113] 

298.15-323.15 

Methanol 

[HMIMJBF4 

DSC 

[157] 

283.15-323.15 

Methanol 

[OMIMJBF4 

Table  IS 

Summary  of  viscosity  measurement  of  alcohol/IL  systems. 

Refrigerant 

Absorbent 

Measuring  method 

Ref. 

Range  of  T/K 

Ethanol 

Methanol 

[DMIMJDMP 

Capillary  viscometer 

[109] 

298.15-323.15 

Methanol 

Methanol 

[OMIMJCI 

[EMIMJDMP 

Capillary  microviscosimeter 

[153] 

298.15-328.15 

[BMIMJDMP 

[EMIMJDMP 

[BMIMJDMP 

Micro  viscometer 

[154] 

293.15-333.15 

Ethanol 

Ethanol 

[EMIM]Tf2N 

Viscosity  meter 

[156] 

293.15-323.15 

■  ,  [BMIM]BF4;  □,  [BMIM]TfO;  A,  [BMIM]PF6;  A,  [EMIM]Tf2N 


Fig.  7.  Schematic  of  the  single-effect  absorption  cooling  cycle. 


The  single-effect  absorption  cooling  cycle  consists  of  an 
evaporator  E,  a  condenser  C,  a  generator  G,  an  absorber  A,  and  a 
solution  heat  exchanger  SHX,  whose  schematic  is  shown  in  Fig.  7. 
1-2-5-3-1  represents  the  refrigerant  circulation  loop,  and 
6-2-5-4-6  is  the  solution  circulation  loop.  In  the  cycle 
calculation,  mass  transfer,  heat  transfer,  and  resistance  loss  in 
the  solution  flow  process  are  neglected.  Therefore,  both  positive 
and  reverse  cycles  can  be  considered  as  ideal  cycles,  with  the 
highest  efficiency  of  heat  transfer  and  cooling  properties. 


In  the  refrigerant  circulation  loop,  the  generator  is  heated  by  Qg 
and  separates  high-pressure  refrigerant  vapor  5.  The  vapor  refrig¬ 
erant  then  enters  the  condenser,  liquefies  into  liquid  refrigerant  3, 
and  releases  heat  Qc.  The  high-pressure  liquid  refrigerant  passes 
through  an  expansion  valve  and  goes  into  the  evaporator.  In  the 
evaporator,  the  liquid  refrigerant  is  vaporized  and  the  cooling  load 
Qe  is  produced. 

In  the  solution  circulation  loop,  the  low-pressure  vapor  refrig¬ 
erant  1  from  the  evaporator  is  absorbed  in  the  absorber  by  the  rich 
solution  8  concentrated  in  the  generator  and  then  releases  heat  Qa. 
The  pressure  of  the  poor  solution  2  leaving  the  absorber  is 
enhanced  after  pumping,  and  then  the  poor  solution  is  supplied 
to  the  generator.  In  the  generator,  the  poor  solution  is  heated  up  to 
boiling  by  the  low-grade  heat  source,  and  the  refrigerant  is 
vaporized  and  separated  from  the  solution.  The  rich  solution 
4  desorbing  the  refrigerant  reduces  the  pressure  through  the 
expansion  valve  and  flows  back  to  the  absorber. 

To  simplify  the  computer  simulation  of  a  cycle,  following 
assumptions  are  made.  The  pressures  of  the  condenser  and  the 
generator  are  equal  (pc=pG),  and  the  pressures  of  the  absorber 
and  evaporator  are  equal  (pa=Pe)-  The  connecting  pipeline  has  no 
pressure  drop.  Compared  with  other  forms  of  thermal  power,  the 
pumping  work  is  so  small  that  is  usually  ignored.  The  refrigerant 
expansion  process  from  the  condenser  to  the  evaporator  is 
isenthalpic,  and  the  temperature  of  evaporator  outlet  is  the  dew 
point  of  the  pure  refrigerant.  The  flow  rate  of  vapor  refrigerant  is 
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set  to  1  kg  s  and  the  vapor  of  absorbent  is  ignored,  etc. 
[158,159], 

According  to  the  above  assumptions,  if  the  refrigerant  vapor 
generated  in  the  evaporator  is  1  kgs-1,  the  heat  added  to  the 
evaporator  is  (/i3+Qe),  and  the  heat  leaving  the  evaporator  is  h2. 
Based  on  the  energy  balance,  the  cooling  capacity  Qe  of  single¬ 
effect  absorption  cooling  cycle  is 

Qe  =  h2  —  h\  (13) 

where  is  the  enthalpy  of  stream  1  in  Fig.  7,  h2  is  the  enthalpy  of 
the  saturated  refrigerant  vapor  at  evaporation  pressure. 

In  the  generator,  when  per  unit  mass  (e.g.,  1  kg)  of  refrigerant 
vapor  is  produced,  /kg  poor  solution  is  needed  to  flow  into  the 
generator.  The  value  of/ is  called  the  circulation  ratio.  Conversely, 
if  the  poor  solution  flowing  into  the  generator  is  / kg,  1  kg 
refrigerant  vapor  is  produced  in  the  generator.  Therefore,  the  rich 
solution  leaving  the  generator  is  (/- 1 )  kg.  The  pure  absorbent 
almost  has  no  vapor  pressure,  i.e.,  the  pure  absorbent  in  the 
stream  inflow  and  outflow  of  the  evaporator  remains  the  same 
based  on  the  mass  balance, 

/wp  =  (/-l)wr  (14) 

or 

/  —  wr/(Wj  Wp)  (15) 

where  wp  is  the  mass  fraction  of  the  absorbent  in  poor  solution, 
wr  is  the  mass  fraction  of  the  absorbent  in  rich  solution. 

Investigating  the  heat  into  and  out  of  the  generator,  the  heat 
added  to  the  generator  is  Qc+fh3,  and  the  heat  leaving  the 
generator  is  (/-l)/i4+fi5,  thus 

Q.c  =  (f~1)h4+h5-fh3  (16) 

The  performance  of  absorption  refrigeration  cycle  is  character¬ 
ized  by  the  coefficient  of  performance  defined  as 

COP  =  Qe/Qg  =  (h2  —  hj  )/[/i5  +f(h4  —  h3)  —  h4]  (17) 

where  h,  represents  the  enthalpy  of  stream  i  (kj  kg  ').  The 
enthalpy  of  each  state  point  can  be  calculated  on  the  basis  of 
thermodynamic  properties  of  refrigerant  or  IL  working  pair. 

The  above  analysis  reveals  that  COP  determination  is  based  on 
the  enthalpy  of  stream.  In  principle,  the  enthalpy  of  stream  is  the 
state  function  on  the  temperature,  pressure,  and  composition  of  a 
given  stream.  In  fact,  Fig.  7  is  also  a  phase  diagram  of  the  binary 
system.  Using  Fig.  7,  the  calculation  method  can  be  explained  by 
the  fact  that  the  state  of  each  point  in  the  cycle  needs  to  be 
determined,  and  then  the  enthalpy  can  be  calculated  based  on  the 
state.  For  example,  based  on  the  heat  source  temperature  Tc,  the 
pressure  and  composition  of  stream  in  the  generator  can  be 
determined.  Generally,  the  enthalpy  of  a  multi-component  fluid  is 


h  =  £ Xjhj+Amixh 

or 

(18) 

rT 

h  =  2>j  /  CPjdT+Amixh 

JTq 

(19) 

where  CPJ  is  the  heat  capacity  of  species  j  of  stream,  Amixh  is  the 
mixing  enthalpy  of  the  system,  which  is  usually  ignored.  T0  is  the 
reference  temperature,  usually  treated  as  298.15  K. 


5.2.  Performance  assessment  of  absorption  cycles 

Table  19  lists  the  published  COP  of  refrigerant/IL  systems  for 
single-effect  absorption  cooling  cycle.  Deviations  among  the 
calculated  results  exist  because  different  calculations  are  based 
on  different  property  data,  thermodynamic  models,  and  operating 
conditions,  e.g.,  heat  source  temperature. 


Table  19 

COP  of  refrigerant/IL  systems  for  single-effect  absorption  cooling  cycle. 


System 

Ref. 

TeI'C 

UC 

TdC 

Tc/C 

/ 

COP 

H20/[BMIM]BF4 

[49] 

10 

30 

40 

100 

13.0 

0.544 

H20/[EMIM]BF4 

18.2 

0.525 

H20/[EMIM]EtS04 

13.57 

0.569 

H20/[DMIM]DMP 

5.32 

0.662 

H20/[BMIM)1 

23.7 

0.534 

H20/[BMIM]DBP 

11.17 

0.532 

H20/[EE1M]DEP 

12.38 

0.565 

H20/[EMIM]DEP 

7.75 

0.588 

H20/[EMIM[DMP 

8.66 

0.691 

H20/[DMIM]DMP 

[72] 

10 

30 

40 

80 

8.77 

0.829 

H20/[DMIM]C1 

9.54 

0.827 

CH40/[DMIM]DMP 

[66] 

10 

30 

40 

100 

3.71 

0.87 

TFE/[BMlM]Br 

[81] 

10 

30 

40 

80 

4.39 

0.807 

TFE/[BMIM]BF4 

6.07 

0.780 

TFE/[EMIM]BF4 

5.25 

0.793 

R22/[BMIM]BFS 

[46] 

10 

30 

40 

100 

5.12 

0.319 

R32/[BMIM]BF6 

7.35 

0.385 

R32/[BM1M]BF4 

6.41 

0.330 

R134/[BM1M]BF6 

4.38 

0.348 

R134a/[BM1M]BF6 

10.66 

0.254 

R152a/[BM1M]BF6 

13.27 

0.300 

R125/[BMIM]BF6 

16.49 

0.128 

NH3/IBM1M]BF4 

[23] 

10 

30 

40 

100 

12.98 

0.557 

NH3/[BMIM]BFs 

17.27 

0.575 

NH3/[EMIM]Tf2N 

24.57 

0.589 

NH3/[HMIM]C1 

14.26 

0.525 

NH3/[EMIM]Ac 

[50] 

10 

30 

40 

100 

12.55 

0.573 

NH3/[EMIM]EtS04 

17.55 

0.485 

NH3/[EMIM]SCN 

12.42 

0.557 

To  study  the  feasibility  of  NH3/IL  systems  in  industrial  applica¬ 
tions,  Yokozeki  et  al.  [50]  calculated  the  performance  of  absorption 
cycle  using  NH3/IL  working  pairs.  In  the  studied  imidazolium  ILs, 
the  cycle  performance  of  [EMIM]Tf2N  as  absorbent  is  the  best,  but 
the  COP  is  slightly  lower  than  that  of  the  NH3/H20  system.  Given 
that  the  distillation  separation  of  the  NH3/[EMIM]Tf2N  system 
does  not  consume  too  much  energy,  the  application  of  NH3/ 
[EMIM]Tf2N  system  may  be  feasible. 

Results  in  Table  19  show  that  the  COP  of  H20/IL  systems  ranges 
within  0.52-0.83,  which  do  not  exceed  the  traditional  H20/LiBr 
working  pair.  However,  the  COP  of  H20/[DMIM[DMP  and  H20/ 
[DM1M]C1  systems  are  close  and  slightly  lower  than  that  of  H20/ 
LiBr.  Particularly,  the  H20/[DM1M]DMP  system  overcomes  the 
crystallization  and  corrosion  problems  of  traditional  working  pair 
to  some  extent,  and  broadens  the  operating  range.  For  HFC/1L 
systems,  the  COP  ranges  within  0.12-0.39,  which  are  still  relatively 
low  but  may  be  used  in  absorption/compression  hybrid  cycle.  For 
NH3/IL  systems,  the  COP  ranges  within  0.48-0.59,  which  are  lower 
than  that  of  the  traditional  NH3/H20  working  pair,  but  the 
difficulty  of  separation  is  relieved. 

Qg  =/(-p(Tg  —  Tjn)+AevahTc  +Amixh  (20) 

where  qG  is  the  heat  load  of  unit  mass  refrigerant  produced  in  the 
generator,  which  consists  of  three  parts.  fCp(Tc-Tin)  expresses  the 
heat  needed  to  warm  the  supercooled  solution  from  entering  the 
generator  to  the  generation  temperature,  Aev3hTC  is  the  latent  heat 
of  evaporation  of  the  refrigerant  at  Tc,  and  Amixh  is  the  mixing 
heat  of  liquid  refrigerant  and  absorbent.  A  larger  /  causes  larger 
heat  load  of  evaporator,  thereby  reducing  the  COP  of  the  cycle. 

Fig.  8  shows  that  with  the  increase  of  generation  temperature, 
the  COP  of  various  working  pairs  sharply  increases  initially,  then 
becomes  stable,  and  slightly  declines  finally.  /  is  infinite  when  the 
generation  temperature  reaches  the  minimum,  the  required  gen¬ 
eration  heat  is  infinite  based  on  Eq.  (20),  and  the  COP  of  the  cycle  is 
zero.  With  the  increase  of  generation  temperature /decreases,  COP 
sharply  increases  and  then  smoothens.  With  further  rise  of 
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Fig.  8.  COP  of  seven  working  pairs  versus  TG  for  single-effect  absorption 
cooling  cycle. 


TJ°C 

Fig.  10.  Effect  of  TE  on  //E. 


Fig.  9.  Effect  of  TE  on  COP. 


generation  temperature,  the  temperature  difference  ( TG-Tin ) 
increases,  the  generated  heat  increases,  and  then  the  COP  of  the 
cycle  slightly  decreases.  At  the  stationary  stage,  the  working  pairs 
consisting  of  TFE  and  three  different  ILs,  whose  COP  ranges 
between  those  of  H20/LiBr  and  NH3/H20  working  pairs,  have 
wider  operating  range.  Among  these  working  pairs,  the  COP  of 
the  TFE/[BMIM]Br  system  is  the  highest.  At  the  stationary  stage, 
the  COP  of  H20/[DM1M]DMP  and  H20/[DMIM]C1  are  close  but 
slightly  lower  than  that  of  the  H20/LiBr  working  pair.  Given  that 
the  H20/[DMIM]DMP  system  has  no  crystallization  limits  and  has 
the  widest  operating  range,  the  operational  safety  greatly 
improves.  The  operating  range  of  the  H20/LiBr  system  comes 
second,  and  the  operating  range  of  the  H20/[DMIM]C1  system  is 
the  narrowest.  As  an  absorbent,  [DMIM]DMP  has  a  higher  solubi¬ 
lity  with  the  refrigerant,  e.g.,  H20,  and  a  higher  boiling  point  than 
refrigerant.  Simulation  results  also  show  that  the  H20/[DMIIVI] 
DMP  system  is  expected  to  become  an  absorption  cycle  alternative 
working  pair. 

For  two  preferred  suitable  H20/IL  systems  H20/[DMIM]C1  and 
H20/[DM1M]DMP,  the  effects  of  evaporation  temperature  on  the 
system  COP,  system  exergy  efficiency  (t/E),  and  /  are  analyzed  and 
compared  with  the  traditional  working  pair  H20/LiBr.  Fig.  9  shows 
the  effect  of  the  evaporation  temperature  on  the  COP,  which  is 
found  to  follow  the  order  H20/LiBr  >  H20/[DMIM]DMP  >  H20/ 
[DM1M]C1,  and  are  all  higher  than  0.78.  With  the  increase  of 
evaporation  temperature,  the  COP  of  the  system  increases,  this  is 
because  that  when  the  evaporation  temperature  is  low,  the 
corresponding  evaporation  pressure  is  low,  /  is  high,  and  COP 
is  low. 

Fig.  10  shows  the  effect  of  the  evaporation  temperature  on  t]E. 
The  effect  of  TE  on  COP  is  opposite  to  that  of  r/E.  COP  increases  with 
the  increase  of  TE,  whereas  rjE  decreases  with  the  increase  of  TE. 


This  is  mainly  because  COP  considers  only  the  amount  of  cooling 
capacity,  while  the  ?/E  considers  not  only  the  amount  of  cooling 
capacity  but  also  the  grade  of  cooling  capacity.  In  other  words,  the 
cooling  load  and  cooling  temperature  are  all  important  for  the 
absorption  cooling  cycle,  a  lower  cooling  temperature  results  in  a 
higher  grade  of  cold,  so  rjE  decreases  with  the  increase  of  TE.  Fig.  10 
shows  that  for  the  three  studied  working  pairs,  i]E  follows  the 
order  H20/LiBr  >  H20/[DMIM]DMP  >  H20/[DMIM]C1. 

Fig.  11  shows  the  effect  of  the  evaporation  temperature  on  /, 
/  values  of  H20/[DMIM]C1  and  H20/[DMIM]DMP  systems  are 
higher  than  that  of  the  H20/LiBr  system,  and  /  decreases  with 
the  increase  of  evaporation  temperature.  /  also  markedly  affects 
COP,  increasing  COP  needs  /  to  be  reduced.  To  reduce  f  the 
concentration  difference  (wr-wp)  should  be  increased  and  the 
concentration  of  the  rich  solution  should  be  reduced. 

As  summarized  above,  the  mixture  systems  of  refrigerant  and 
IL  are  novel,  environmental  friendly  absorption  cycle  working 
pairs.  Thus,  these  systems  are  attracting  increasing  attentions, 
and  further  research  and  development  are  necessary. 


6.  Conclusions 

The  research  frontiers  of  IL  working  pairs  for  absorption  cycle 
in  recent  years  were  surveyed,  to  make  the  relevant  work  more 
rational  and  efficient.  Based  on  macroscopic  properties  and  inter- 
molecular  interactions,  selection  methods  of  IL  as  absorbent  were 
introduced.  In  particular,  based  on  the  author's  previous  studies, 
from  the  perspective  of  action  mechanism  between  absorbent  and 
refrigerant,  evaluation  methods  of  combining  the  UNIFAC  model 
with  extreme  GE  criterion  were  proposed,  and  the  actual  research 
achievements  of  the  method  were  also  presented.  The  research 
progresses  of  thermophysical  properties  of  IL  working  pairs  were 
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summarized,  such  as  the  test  method,  popular  system,  modeling  of 
vapor  pressure,  solubility,  heat  capacity,  and  density.  Based  on 
thermophysical  property  investigations  of  H20,  NH3,  alcohol,  and 
HFC  systems,  as  well  as  researches  of  ternary  systems  with  IL 
added  to  the  traditional  working  pairs,  several  potential  binary 
and  ternary  systems  were  put  forward.  Finally,  working  perfor¬ 
mances  of  absorption  cycle  with  IL  working  pairs  were  calculated 
and  assessed  to  further  validate  the  development  potential  of  H20, 
NH3,  alcohol,  and  HFC  working  pairs. 


Nomenclature 

At  and  parameters  of  Eq.  (10)  (dimensionless) 

An,  Bn  and  Cn  parameters  of  Eq.  (12)  (dimensionless) 
a,  b,  and  c  parameters  of  Eq.  (3)  (dimensionless) 

dj  and  fa,  parameters  of  Eq.  (11)  (dimensionless) 

CFC  chlorofluorocarbon  (dimensionless) 

COP  coefficient  of  performance  (dimensionless) 

Cp  the  heat  capacity  (kj  kg-1  K_1) 

DBP  dibutylphosphate  (dimensionless) 

DMEU  N,N'-dimethylethyleneurea  (dimensionless) 

DMPU  N,N'-dimethylpropyleneura  (dimensionless) 

DSC  differential  scanning  calorimetry  (dimensionless) 

E181  tetraethylene  glycol  dimethyl  ether  (dimensionless) 

EOS  equation  of  state  (dimensionless) 

/  circulation  ratio  (dimensionless) 

GE  excess  Gibbs  function  (kj  kg-1) 

HC  hydrocarbon  (dimensionless) 

HFC  hydrofluorocarbon  (dimensionless) 

HFIP  hexafluoroisopropanol  (dimensionless) 

hi  enthalpy  of  stream  i  in  Fig.  7,  i=  1,2 . 5  (kj  kg-1) 

IL  ionic  liquid  (dimensionless) 

k i  Henry's  constant  of  the  refrigerant  species  (kPa) 

NMP  N-methyl-2-pyrrolidone  (dimensionless) 

p  system  pressure  (kPa) 

Pi  saturated  vapor  pressure  of  the  refrigerant  species  (kPa) 

Q  heat  load  (kW) 

q  heat  load  of  unit  mass  refrigerant  (kW  kg-1) 

R&D  research  and  development  (dimensionless) 

R114  l,2-dichloro-l,l,2,2-tetrafluoroethane  (dimensionless) 

R125  pentafluoroethane  (dimensionless) 

R134  1,1,2,2-tetrafluoroethane  (dimensionless) 

R134a  1,1,1,2-tetrafluoroethane  (dimensionless) 

R14  tetrafluoromethane  (dimensionless) 

R143a  1,1,1 -trifluoroethane  (dimensionless) 

R152a  1,1-difluoroethane  (dimensionless) 

R161  fluoroethane  (dimensionless) 

R23  trifluoromethane  (dimensionless) 

R32  difluoromethane  (dimensionless) 

R41  fluoromethane  (dimensionless) 

SHX  solution  heat  exchanger  (dimensionless) 

T  system  temperature  (K) 

T0  reference  temperature  (K) 

TFE  2,2,2-trifluoroethanol  (dimensionless) 

Tin  temperature  of  supercooled  solution  entering  the 

generator  (K) 

VLE  vapor-liquid  equilibrium  (dimensionless) 

w  mass  fraction  of  the  absorbate  species  (dimensionless) 

x,  liquid  phase  mole  fraction  of  species  i  (dimensionless), 
i=l,  2,  1  presents  refrigerant,  2  presents  absorbant 

y, -  vapor  phase  mole  fraction  of  species  i  (dimensionless), 

i=l,  2,  1  presents  refrigerant,  2  presents  absorbant 
[BMIM]BF4  l-butyl-3-methylimidazolium  tetrafluoroborate 
(dimensionless) 


[BMlM]Br  l-butyl-3-methylimidazolium  bromide 
(dimensionless) 

[BM1M]C1  l-butyl-3-methylimidazolium  chloride 
(dimensionless) 

[BM1M]DBP  l-butyl-3-methylimidazolium  dibutylphosphate 
(dimensionless) 

[BMlM]MeS04  l-butyl-3-methylimidazolium  methylsulfate 
(dimensionless) 

[BM1M]PF6  l-butyl-3-methylimidazolium  hexafluorophosphate 
(dimensionless) 

[BMlM]Tf2N  l-butyl-3-methylimidazolium  bis(trifluoromethyl- 
sulfonyl)imide  (dimensionless) 

[BMlM]TfO  l-butyl-3-methylimidazolium  trifluoromethanesulfo- 
nate  (dimensionless) 

[C10MlM]Br  l-decyl-3-methylimidazolium  bromide 
(dimensionless) 

[C10M1M]C1  l-decyl-3-methylimidazolium  chloride 
(dimensionless) 

[CnMIM]BF4  l-all<yl-3-methylimidazolium  tetrafluoroborate, 
n  =  l,  2,...  (dimensionless) 

[CnMlM]Cl  l-all<yl-3-methylimidazolium  chloride,  n  =  l,  2,... 
(dimensionless) 

[CnMlM]DMP  l-alkyl-3-methylimidazolium  dimethylphosphate, 
n  =  l,  2,...  (dimensionless) 

[CnMlM]EtS04  l-alkyl-3-methylimidazolium  ethylsulfate,  n=l, 
2,...  (dimensionless) 

[CnMlM]MeS04  l-alkyl-3-methylimidazolium  methylsulfate, 
n=l,  2,...  (dimensionless) 

[CnMIM]PF6  l-alkyl-3-methylimidazolium  hexafluorophosphate, 
n=l,  2,...  (dimensionless) 

[CnMIM]Tf2N  l-alkyl-3-methylimidazolium  bis(trifluoromethyl- 
sulfonyl)imide,  n  =  1,  2,...  (dimensionless) 

[CnMIM]TfO  l-alkyl-3-methylimidazolium  trifluoromethanesul 
fonate,  n  =  1,  2,...  (dimensionless) 

[DMEA]Ac  N,N-dimethylethanolammonium  acetate 
(dimensionless) 

[DM1M]BF4  l-methyl-3-methylimidazolium  tetrafluoroborate 
(dimensionless) 

[DM1M]C1  1,3-dimethylimidazolium  chloride 
(dimensionless) 

[DM1M]DMP  1,3-dimethylimidazolium  dimethylphosphate 
(dimensionless) 

[EEIM]DEP  l-ethyl-3-ethyl-imidazolium  diethylphosphate 
(dimensionless) 

[EMIM]BEI  l-ethyl-3-methylimidazolium  bis(pentafluoroethyl 
sulfonyl)imide  (dimensionless) 

[EMIM]BF4  l-ethyl-3-methylimidazolium  tetrafluoroborate 
(dimensionless) 

[EMlM]Br  l-ethyl-3-methylimidazolium  bromide 
(dimensionless) 

[EMIM]Cll-ethyl-3-methylimidazolium  chloride 
(dimensionless) 

[EMIM]DEP  l-ethyl-3-methylimidazolium  diethylphosphate 
(dimensionless) 

[EM1M]DMP  l-ethyl-3-methylimidazolium  dimethylphosphate 
(dimensionless) 

[EMIM]EtS04  l-ethyl-3-methylimidazolium  ethylsulfate 
(dimensionless) 

[EMIM]Tf2N  l-ethyl-3-methylimidazolium  bis(trifluoromethyl- 
sulfonyl)imide  (dimensionless) 

[EMIM]TfO  l-ethyl-3-methylimidazolium  trifluoromethanesulfo- 
nate  (dimensionless) 

[HMIM]BF4  l-hexyl-3-methylimidazolium  tetrafluoroborate 
(dimensionless) 

[HMIM]C1  l-hexyl-3-methylimidazolium  chloride 
(dimensionless) 
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[HMIM]PF6  l-hexyl-3-methylimidazolium  hexafluorophosphate 
(dimensionless) 

[HMIM]Tf2N  l-hexyl-3-methylimidazolium  bis(trifluoromethyl- 
sulfonyl)imide 
(dimensionless) 

[HMIM]TfO  l-hexyl-3-methylimidazolium  trifluoromethanesul- 
fonate  (dimensionless) 

[HydeMIM]BF4  l-(2-hydroxyethyl)-3-methylimidazolium  tetra 
fluoroborate  (dimensionless) 

[Nih(2oh)]C1  choline  chloride  (dimensionless) 

[P4444]TMPP  tetrabutylphosphonium  bis(2,4,4-trimethylpentyl) 
phosphinate  (dimensionless) 

[P8111]TMPP  trimethyloctylphosphonium  bis(2,4,4-trimethyl 
pentyl)phosphinate  (dimensionless) 

Creek  Letters 


Yi 

activity  coefficient  of  species  i  (dimensionless), 

1  presents  refrigerant,  2  presents  absorbant 

i=  1,  2, 

•^eva^rc 

latent  heat  of  evaporation  of  the  refrigerant 
(kJkg  ’) 

at  Tc 

Amixh 

mixing  enthalpy  of  the  system  (kj  kg  ') 

IE 

system  exergy  efficiency  (dimensionless) 

P 

density  (g  cm~3) 

¥ri 

absorption  potential  of  the  refrigerant 
(dimensionless) 

species 

Superscripts 

C 

combinatorial  term  (dimensionless) 

R 

residual  term  (dimensionless) 

RL 

Raoult’s  law  (dimensionless) 

oo 

infinite  dilution  state  (dimensionless) 

Subscripts 

A 

absorber  (dimensionless) 

C 

condenser  (dimensionless) 

E 

evaporator  (dimensionless) 

G 

generator  (dimensionless) 

max 

maximum  value  (dimensionless) 

P 

poor  solution  (dimensionless) 

r 

rich  solution  (dimensionless) 
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